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I. INTRODUCTION 


Reference 1 presents the electromagnetic theory for a thin-wire 
structure in a homogeneous conducting medium, and this report presents 
the corresponding computer program. The program performs a frequency- 
domain analysis of thin-wire antennas and scatterers. The wire con- 
figuration is a generalized polygon assenbled from straight wire seg- 
ments. The program has been tested extensively with simple structures 
(linear dipoles, V dipoles, coupled dipoles, square. loops, octagonal 
loops, multi turn loops and coupled loops) and complicated configura- 
tions incl uding wi re-gri d models of plates, spheres, cones, aircraft 
and ships. Although the air-earth or air-water interface is not con- 
sidered, the program is applicable in many problems involving buried 
or submerged antennas or targets. It is useful in locating the poles 
of the admittance or scattering function for wire structures in the 
complex frequency domain. 

A piecewi se-sinusoi dal expansion is used for the current distri- 
bution. The matrix equation ZI = V is generated by enforcing reaction 
tests with a set of sinusoidal dipoles located in the interior region 
of the wire. Since the test dipoles have the same current distribution 
as the expansion modes, this may be regarded as an application of 
Galerkin's method. Rumsey's reaction concept was most helpful in this 
development, and therefore the formulation is known as the "sinusoidal 
reaction technique." 

The current is assumed to vanish at the endpoints (if any) of the 
wire, and Kirchhoff's current law is enforced everywhere on the 
structure. The input data specify the frequency, wire radius, wire 
conductivity, the parameters of the exterior medium, coordinates of 
points to describe the shape and size of the wire configuration, and a 
list of the wire segments. If some or all of the wire segments are 
insulated, the radius and permittivity of the insulating sleeve are 
indi cated. 

Coordinates are required for wire endpoints, corners, junctions 
and terminals. For accuracy, the longest wire segment should not 
greatly exceed one-quarter wavelength. Longer segments should be 
subdivided by defining additional current-sampling. points. The pro- 
gram automatically defines a set of N sinusoidal dipole modes on the 
wire structure and computes the mutual impedance matrix for these 
modes. The elements in the matrix are generated by numerical inte- 
gration when appropriate, or from closed-form expressions in terms of 
exponential integrals. The computer program uses certain approxima- 
tions which yield a symmetric matrix even when the wire structure has 
finite conductivity, lumped loads and insulating sleeves. 

In antenna problems, the output data Includes the current distri- 
bution, impedance, radiation efficiency, gain, patterns and near-zone 
field. In bistatic scattering problems, the output includes the echo 


area and the complex elements of the polarization scattering matrix. 

In backscatter situations, the output includes also the absorption, 
scattering and extinction cross sections. 

If the wire has finite conductivity or dielectric sleeves, it is 
assumed that the frequency is real. This restriction can readily be 
removed if the user will specify the surface impedance of the wire and 
the complex permittivities of the dielectric sleeves and the ambient 
medium appropriate for complex frequencies. 

The user may make a tradeoff between accuracy and computation 
costs by specifying the input variable INT. A large value Increases 
the accuracy and the cost. For most problems, the reconmended value 
is INT = 4. 

The program was run on an IBM 370/165 computer to determine the 
broadside backscatter for a wire-grid square plate with edge length L. 
With a five-by-five grid, there are 60 segments, 36 points and 84 
simultaneous equations. With INT = 4, calculations were made for L/x 
= 0.3, 0.4, 0.5, 0.6 and 0.7. The execution time was 100 seconds. 

This averages to 20 seconds for each value of L/\. The wire structure 
was perfectly conducting, uninsulated and located In free space. No 
advantage was taken of the target symmetries. 

The next section presents the thin-wire computer program, instruc 
tions for the user, typical input and output data and tables of the 
mutual impedance of sinusoidal dipoles. Appendicles list the computer 
subroutines and explain their functions. 

II. THE TNIN-WIRE COMPUTER PROGRAM 

Fig. 1 is a Fortran listing of the thin-wire computer program. 
Near the beginning of this program, the DIMENSION statements reserve 
storage for a wire structure with up to 50 segments, 55 points and 60 
dipole modes. Quantities with the same or related dimensions are 
grouped together on the same line or consecutive lines. 

NM denotes the actual number of monopoles (segments), INM is the 
corresponding dimension, and the dimension for CG, VG and ZLD Is twice 
INM. The second subscript for MD always has a dimension of 4. 

N denotes the nunber of simultaneous linear equations and ICJ is 
the corresponding dimension. The dimension for C is (ICJ*ICJ + ICJ)/2 

The DO LOOP ending at statement 15 sets ISC(J) = 0 for all the 
segments. This indicates that the wires are bare or uninsulated. If 
some or all of the segments are insulated, the user may set ISC(J) = 1 
for the appropriate segment numbers J. 
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COMPLEX EP2,EP3,ETA,GAM,Y11,Z11,ZS OOOl 

COMPLEX EPPS,EPTS,ETPS,ETTS,EX,EY,EZ 0002 

COMPLEX C(1830),CJ(60),EP(60),ET(60),EPP(60),ETT(60) 0003 

DIMENSION 11(60), 12(60), I 3 ( 60 ) » JA( 60 ) » JB( 60) 0004 

COMPLEX CGD (50) ,SGD (50 ) ,CG! 100) ,VG( 100) ,ZLD( 100) 0005 

DIMENSION D(50),IA(50),IB(50),ISC(50),MD(50,4),ND(50) 0006 

DIMENSION X(55 ) » Y ( 55 ) , Z ( 55 ) 0007 

DATA PI , TP/3. 14159 ,6.28318/ 0008 

DATA EO »U0/ 8 . 854E -1 2 , 1 .2566E-6/ 0009 

2 FORMAT! IX ,8F 15 .7 ) 0010 

3 FORMAT! 1X»4F15.7/) 0011 

4 FORMAT! IX, 115, 8F14. 6) 0012 

5 FORMAT ( 1 HO ) 0013 

6 FORMAT! 1X,6F15.7/) 0014 

7 FO RM AT ( 8F 10 .5 ) 0015 

8 FORMAT! IX, 114, 1315 ) 0016 

9 F0RMATI3X, 'MAX = ',I5,3X,'MIN = ',I5,3X,'N « ',151 0017 

I C J = 60 0018 

INM =5 0 0019 

DO 15 J = 1 , INM OOZO 

15 ISC(J)=0 0021 

RE AD ( 5 ,7 ) BM ,ER2 ,S 1 G2 ,TD2 0022 

WRITE (6,2)BM,ER2,SIG2,TD2 0023 

RE AD ( 5 ,7)AM,CMM,ER3,SIG3,TD3 0024 

WRITE (6,2) AM, CMM,ER3,SIG3,TD3 0025 

RE AD ( 5 , 8 ) I B 1 SC , I G A IN , INE AR , I SC AT , IW R ,N GEN , NM, NP 0026 

WRITE (6,8) IBISC,IGAIN, INE AR , I SC AT , IWR ,NGEN, NM , NP 0027 

RE AD (5 ,7)FMC»PHA,THA»PHI ,THI ,PHS,THS 002 8 

WRITE ! 6,2 )FMC , PHA , T H A , PH I ,THI ,PHS, THS 00 29 

DO 22 J = 1 ,NM 0030 

REA0(5,8)IA(J),IB!J) 0031 

22 WRITE ( 6,8 ) J , I A! J ) , IB! J ) 0032 

DO 40 1=1, NP 0033 

RE AD ( 5 , 7 ) X ( I ) , Y ( 1 ) , Z ( I ) 0034 

40 WRITE (6,4) I ,X( I ) ,Y( I ) ,Z( I ) 0035 

READ(5,7)XP,YP»ZP 0036 

FHZ=FMC*1.E6 0037 

OME GA = TP*F HZ 0038 

IF (SIG2.LT.0.)EP2=ER2*E0*CMPLX(1.,-TD2) 0039 

IF (TD2.LT.0. )EP2=CMPLX(ER2*E0»-S1G2/0MEGA> 0040 

IF (S1G3.LT.0.)EP3=ER3*E0*CMPLX!1.,-TD3) 0041 

IF (TD3.LT.0. )EP3=CM PLX(ER3*E0, -SI G3 /OMEGA) 0042 

ETA=CSQRT!U0/EP3) 0043 

GAM =UM EGA *C SORT! -U0*EP3) 0044 

CALL SORT! IA, IB, I 1,12, 13, JA,JB,MD ,ND ,NM ,N P , N,MAX , MI N , I C J, INM) 004 5 

W R I TE ( 6 ,5 ) 0046 

WRITE (6,9)MAX,MIN,N 0047 

W R I TE ( 6 ,5 ) 0048 

IFIMAX.GT.4 .OR. MIN.LT.l .OR. N.GT.1CJ)G0 TO 800 0049 

IN T =4 OO 50 

112=1 005 1 

DO 60 J=1 »NM 0052 

VG ( J ) = ( .0, .0 ) 0053 

ZLD(J) = ( • 0 , . 0 ) 0054 

J J = J+NM 0055 

VG! JJ ) = ( .0 , .0 ) 0056 

60 Z LD ( J J ) = ( .0 , .0 ) 0057 

IF (NGEN.GT.O )VG!NGEN )=( 1. ,0. ) 0058 

CALL SGANT(IA,IB,INM,INT,ISC,I1,I2,I3,JA,JB,MD,N,ND,NM,NP 0059 

2*AM,BM,C,CGD ,CMM ,D ,E P2 ,E P3 ,E T A ,F HZ ,G AM , SGD, X ,Y , Z , Z L D, Z S) 0060 

IF (N.LE.O)GO TO 800 006 1 

IF (NGEN.LE .0)G0 TO 400 °0*>2 


Fig. la. The thin-wire computer program. 
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CALL GANT1 ( I A , I B , INM , IWR , 1 1 , 1 2 , 1 3 , 1 12 , JA, JB> MD,N, ND, NM, AM 0063 

2. C,CJ,CG,CMM ,D ,EFF ,GAM ,GG,CGD , SGO ,VG ,Y 11 ,Z 1 1 » ZLD, ZS ) 0064 

WRITE (6,3)EFF ,GG,Z11 0065 

200 IF ( INEAR. LE .0)G0 TO 300 0066 

CALL GNFLO (I A, IB ,INM ,1 1,1 2,13 ,MD,N ,ND , NM , AM, CGO , SGO, ETA, GAM 0067 

2,CJ,D,X,Y,Z,XP,YP,ZP,EX,EY,EZ) 0068 

WRITE (6,3)XP,YP,ZP 0069 

WRITE (6,6JEX,EY,EZ 0070 

300 IF I1GAIN.LE .0)G0 TO 400 0071 

INC=0 0072 

PH=PHA 0073 

TH=THA 0074 

CALL GFFLO ( I A , IB ,INC ,INM , IWR,I1,I 2,13,1 12,MD,N,ND,NM,AM 0075 

2, ACSP,ACST,C,CGD ,CG,CJ,CMM ,D ,E CSP ,E C S T ,E P , ET, EPP , E T T, EPPS, EPTS 0076 

3, ETPS,ETTS ,GG ,GPP , GTT , PH , SGO , SCSP » SCST, SPPM , SP TM f STPM , S TTM, TH 0077 

4, X,Y,Z,ZL0,ZS,ETA,GAM) 0078 

WRITE 16, 3)PH,TH, GPP, GTT 0079 

400 IF (ISCAT.LE.O)GO TO 600 0080 

INC= 1 0081 

PH= PH I 0082 

TH=TH I 0083 

CALL GFF LD( I A, IB, INC, INM ,IWR,I 1,12,1 3,1 12 , MD, N, ND, NM, AM 0084 

2»ACSP,ACST,C,CGD,CG,CJ,CMM ,D ,ECSP ,E CST,E P , ET , EPP , ETT, EPPS , EPTS 0085 

3, ETPS,ETTS,GG,GPP ,GTT ,PH,SGD ,S CSP , SCST, SPPM, SPTM, STPM, STTM, TH 0086 

4, X,Y,Z,ZLD ,ZS , ETA, GAM ) 0087 

WRITE (6,6 ) PH, TH, SPPM , SPTM , STPM , STTM 0088 

WRITE (6,6) ACS P , ACST , ECSP ,E CST , SCSP, SCST 0089 

500 I F ( I8ISC.LE.01GO TO 600 0090 

INC=2 0091 

PH=PHS 0092 

TH=THS 0093 

CALL GFF LD ( I A, I B, INC, INM ,IWR ,1 1, 12, 1 3, 1 1 2 ,MD » N, ND, NM, AM 0094 

2»ACSP,ACST,C,CGD ,CG,CJ ,CMM ,0 ,ECSP ,E CST.E P , ET , EPP , ETT , EPPS , E PTS 0095 

3, ETPS,ETTS,GG,GPP ,GTT ,PH,SGD ,S CSP , SC ST, SPPM, SPTM, STPM, STTM, TH 0096 

4, X,Y,Z,ZLD ,ZS , ETA, GAM) 0097 

WRITE (6, 6)PH,TH, SPPM , SPTM , STPM , STTM 0098 

600 CONTINUE 0099 

800 CALL EXIT 0100 

END 


Fig. lb. The thin-wire computer program. 
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The first READ statement inputs the following parameters for the 
dielectric insulation: 

BM outer radius in meters 

ER2 dielectric constant relative to free space 

SIG2 conductivity in mhos per meter 

TD2 loss tangent 

The program will use SIG2 or TD2 but not both. The user determines 
which one will be used by assigning the other a negative value. For 
an uninsulated wire structure, the program will not use any of the 
data from the first READ statement. 

The second READ statement inputs the following parameters for the 
wire and the exterior medium: 

AM wire radius in meters 

CMM wire conductivity in megamhos per meter 

ER3 dielectric constant relative to free space 

SIG3 conductivity in mhos per meter 

TD3 loss tangent 

The parameters ER3, SIG3 and TD3 are those of the homogeneous ambient 
medium. Again, the program will use SIG3 or TD3 but not both. 

The third READ statement inputs the following data: 

IBISC indicator for bistatic scattering calculations 
I GAIN indicator for antenna gain calculations 
INEAR indicator for near-zone field calculations 

I SCAT indicator for backs catter calculations 
IWR indicator for write out of current distributions 

NGEN indicator for antenna calculations 

NM number of monopoles (segments) 

NP number of points 

For each indicator, a positive value means the calculation or writeout 
is desired while a zero or negative value means it is not desired. 

The fourth READ statement inputs the following data: 

FMC frequency in megahertz 

PHA,THA far- field angle for antenna gain 

PHI ,THI incidence angle for plane-wave scattering 

PHS.THS scattering angle for bistatic scattering 

The above angles are given in degrees, and they denote values of the 
angular coordinates in the spherical system (r,9,<ti) widely used in 
antenna and scattering literature. 
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The fifth READ statement (in the DO LOOP ending with statement 22) 
inputs the endpoints IA(J) and IB(J) of segment J. Thus, IA and IB are 
the Index numbers of the two points which are joined by segment J. 

The sixth READ statement (in the DO LOOP ending with statement 40) 
inputs the coordinates X(I), Y(I) and Z(I) of point I in meters. The 
seventh and last READ statement inputs the coordinates XP, YP and ZP 
(in meters) of the observation point for near-zone field calculations. 

Some of the quantities used in the program are defined as follows: 

FHZ frequency in Hertz 

OMEGA angular frequency 

EP2 complex permittivity of insulation 

EP3 complex permittivity of ambient medium 

ETA intrinsic impedance of ambient medium 

GAM intrinsic propagation constant of ambient medium 

ZS surface impedance of wire 

For an uninsulated wire with perfect conductivity, one may specify 
complex values for ETA and GAM and delete the following input data and 
calculations: BM, ER2, SIG2, TD2, ER3, SIG3, TD3, FMC, FHZ, OEMGA, 

EP2 and EP3. 

After reading the input data, the program calls subroutine SORT. 
This subroutine defines a set of dipole modes on the wire structure. 

N denotes the total number of dipole modes, the number of simultaneous 
linear equations, and the size of the impedance matrix Z. .. Since this 
matrix is symmetric, only the upper-right triangular portion (including 
the entire principal diagonal) is calculated and stored in C(K). SORT 
calculates N, but the user may predict N as follows to reserve adequate 
storage. If m wire segments intersect at a point, this point is de- 
fined as a junction of order m and degree n = m - 1. There will be n 
dipole modes with terminals at this junction. N is determined by 
summing the degrees of all the junctions. For an example, an endpoint 
of a dipole is a junction of order m = 1 and degree n = 0. The vertex 
of a V dipole is a junction of order 2 and degree 1. NP denotes the 
number of points on the wire structure, and each of these points is 
considered to be a junction. 

Mode I is a two-segment V dipole with a sinusoidal current dis- 
tributed over the intersecting segments JA(I) and JB(I). The dipole 
has endpoints 11(1) and 13(1) and terminals at 12(1). The reference 
direction for positive current on dipole I is from II to 12 to 13. 

A wire segment may be shared by as many as four dipole modes, or 
as few as one. In the output of subroutine SORT, ND(J) denotes the 
number of dipoles sharing segment J. The extreme values of ND( J ) are 
MAX and MIN. If MIN is less than one, the wire structure has an un- 
connected segment and the computation is aborted. (An isolated wire 
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must have at least two segments and three points.) If N exceeds ICJ, 
the dimensions are inadequate and the run Is aborted. 

I NT specifies the number of Intervals for calculating the elements 
in the impedance matrix with Simpson's-rule Integration. A large value 
for INT Improves the accuracy at the expense of greater execution time. 

For most problems a suitable combination of speed and accuracy is ob- 
tained with INT = 4. A larger value is recommended if one wire passes 
close to another as in the helix or the multitum loop. If in doubt, 
one may set INT = 0 to choose the rigorous closed-form impedance ex- 
pressions in terms of exponential integrals. 

The DO LOOP ending with statement 60 sets all the lumped load 
impedances and generator voltages to zero. If the wire structure has 
lumped loads, one may insert a READ command after statement 60 to input 
a list of complex load impedances ZLD(J). For a wire antenna with just 
one generator, the program inserts a unit voltage generator wi th VG(NGEN) 
= (l.,0.). If the antenna or array has several generators, one may 
insert a READ command after statement 60 to input a list of complex 
voltages VG(0). 

Generators or lumped loads may be inserted at either end or both 
ends of segment J. First consider a load impedance inserted in the 
middle of segment J. Now slide the load along the segment and let it 
approach the endpoint IA(J). This load is represented by ZLD(J). Next 
insert another load in segment J and slide it to approach the endpoint 
IB(J) . This load is designated ZLD(JJ) where JJ = J + NM. The same 
convention is employed for the voltage generators VG(J) and VG(JO). 

A generator voltage VG(J) is considered positive if it tends to force 
a current flow in the direction from I A ( J) to IB(J). 

Subroutine SGANT calculates the elements of the impedance matrix 
1 - • and stores them in C(K) where K = ( I - 1 ) *N - (1*1 - I)/2 + J. This 
subroutine will set N = 0 and the run will abort if the wire radius is 
zero or negative, the shortest segment length is less than the wire 
diameter, the wire radius is electrically large, or the longest segment 
is too long. 

Subroutine GANT1 considers the thin-wire structure as an antenna 
and solves for the current distribution CG(J), radiation efficiency 
EFF, time-average power input GG and complex power input Yll. In the 
current distribution, CG(J) is the current on segment J as one approaches 
the endpoint IA(J) and CG(JJ) is the current at the other end IB(J). 

The reference direction for positive current is from IA to IB. Thus, 
the conventions are the same for the branch currents CG and the branch 
voltages VG. 

If the antenna has only one voltage generator with VG(NGEN) = 
(l.,0.), then Yll is the antenna admittance and Zll is the impedance. 
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The radiation efficiency EFF is calculated from the time-average 
power Input to the antenna and the time-average power dissipated in 
the wire and the lumped loads. If the antenna Is insulated, the power 
dissipated in the insulation is neglected. If the wire has perfect 
conductivity and the loads are purely reactive, the calculated effi- 
ciency will be 100 per cent. 

The near-field subroutine GNFLD calculates the electric field 
intensity (EX,EY,EZ) at the observation point (XP,YP,ZP). In the 
calling parameters, CJ denotes the current distribution on the wire. 

(The loop currents are stored in CJ(I) and the branch currents in CG(J)). 
Thus, the currents must be calculated before GNFLD is called. Fig. 1 
illustrates the use of GNFLD to calculate the near-zone field In an 
antenna problem. This subroutine can be called again just above 
statement 500 to calculate the near-zone scattered field for a wire 
target. In the calling parameters, CJ is replaced with EP or ET to 
obtain the near-zone field with a phi -polarized or theta-polarized 
incident plane wave. Reference 1 describes the more sophisticated 
techniques required when the observation point is extremely close to 
the wire structure. 

The far-field subroutine GFFLD calculates antenna gain if INC = 0, 
backscattering if INC = 1, and bistatic scattering if INC = 2. If 
INC = 0, PH and TH denote the spherical coordinates <► and 0 of the 
distant observation point and the output from GFFLD is defined as 
follows. EPPS and ETTS denote the phi -polarized and theta-polarized 
components of the electric field intensity. For exanple, 

(1) EPPS = r e yr E a 

9 

where r is the distance from the origin to the observation point. GPP 
and GTT denote the power gains associated with the phi and theta 
polarizations. Appendix 14 defines GPP and GTT more precisely. 

If INC = 1, PH and TH denote the incidence angles 9 - and e.. 

These are also the spherical coordinates of the distant iource. 1 In 
this backscattering situation, the output data from GFFLD are defined 
as follows: 

ACSP,ACST 
ECSP,ECST 
EP.ET 
EPPS 
EPTS 
ETPS 
ETTS 

SCSP,SCST 
SPPM 


absorption cross sections for 9 and e polarizations 

extinction cross sections for 4 and 0 polarizations 

loop currents induced by 9 and e polarized waves 

scattered electric field E a. 

scattered electric field E* e 

scattered electric field Eg. 

scattered electric field E Q ! 

scattering cross sections for 9 and 0 polarizations 
echo area cr. i 
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SPTM echo area (Taq 

STPM echo area 

STTM echo area cr && 

The echo areas are given in square meters. For the doubly-subscripted 
quantities such as E- and tr**, the first and second subscripts specify 
the polarizations of the incident and scattered waves, respectively. 

The complex numbers EPPS, EPTS, ETPS and ETTS are the elements of the 
polarization scattering matrix. 

If INC = 2, PH and TH denote the scattering angles <|> s and 8 $ . 

These are the spherical coordinates of the distant observer. In this 
bi static scattering situation, the only outputs from GFFLD are the 
polarization scattering matrix and the echo areas. 

To obtain antenna patterns, backscattering patterns or bistatic 
patterns, one may insert DO LOOPS in the program to increment the 
angles PH and TH. The DO LOOP will begin just above the call to GFFLD 
and terminate just below this call. To obtain the near-zone field 
distribution along a given probing path, one may insert a DO LOOP 
beginning just above the call to GNFLD and terminating just below this 
call . 


When the calculations have been completed for one problem, one may 
GO TO a point just above CALL GANT1 If only the generator voltages are 
to be changed. One may GO TO a point just below CALL SORT if there is 
a change in the wire radius or conductivity, the insulation, ambient 
medium, frequency, load impedances or the coordinates (X,Y,Z). If 
there is a change In NM, NP, IA or IB, one should GO TO a point above 
CALL SORT. 

Consider an array of three center-fed dipoles, and suppose we 
desire the 3 x 3 admittance matrix for the array. Let each dipole be 
divided into four segments with segments 1 through 4 on dipole 1, 5 
through 8 on dipole 2 and 9 through 12 on dipole 3. The three-port 
admittance matrix can be obtained by inserting a DO LOOP beginning 
just above CALL GANT1 and terminating just below this call. GANT1 will 
be called three times with all the voltages VG set to zero except for 
a single one-volt generator. On the first, second and third calls, 
let NGEN = 3, 7 and 11 to represent a generator at the center of dipole 
1, 2 and 3, respectively. After the first call, set Yll = CG(3), 

Y12 = CG(7) and Y13 = CG(ll). Set Y22 = CG(7) and Y23 = CG(ll) after 
the second call and Y33 = CG(ll) after the third call. 

For extremely small antennas, quasi-static or double-precision 
subroutines are required. 

The wire radius must exceed zero, but there is no difficulty with 
small radii. If the radius exceeds 0.007X, the thin-wire assumptions 
are questionable and the accuracy and convergence deteriorate. The length 
ratio of the longest and shortest segments should not exceed 100. It is 
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assumed that the wire length exceeds the wire diameter by a factor of 
at least 30. We are not aware of any lower limit on the segment length, 
however. 

If a wire is bent sharply to form a small acute angle (less than 
30 degrees), the thin-wire model is questionable. It Is assumed that 
the wire conductivity greatly exceeds the conductivity of the ambient 
medium. For insulated wires, the dielectric layer Is assumed to be 
electrically thin. 

For each thin-wire problem, calculations should be repeated several 
times with the wire divided progressively into shorter segments. There 
is no assurance of accuracy until the output data converge. For a 
moderately thick wire (with radius a = 0.007 \ or larger), the suscept- 
ance may diverge with the delta-gap model. This difficulty may be 
alleviated or eliminated with the magnetic-frill model and the tech- 
niques of Imbrlale and Ingerson [Z ] . 

Tables 1, 2 and 3 list input and output data for three simple 
examples of uninsulated wire structures.. Each table includes a sketch 
of the wire configuration with labels to Indicate the nunfcering system 
for the points and segments. In these examples there are no lumped 
1 oads . 

In the sinusoidal -reaction formulation, a basic function is the 
mutual impedance between two sinusoidal filamentary electric dipoles. 

One dipole is a test source located on the axis of the wire structure, 
and the other is an expansion mode on the wire surface. In view of the 
importance of this mutual impedance, short tables are presented next 
for a few simple cases. The data can be reproduced with the program 
In Fig. 1 with appropriate Input data for uninsulated wires with perfect 
conductivity and no lumped loads in free space. The data were obtained 
with the closed-form expressions (INT = 0) by writing 'out the quantities 
C(K) just below the call to subroutine SGANT. Double precision was 
used for these calculations. 

Table 4 lists the self impedance of a two-segment sinusoidal V 
dipole with radius a = 0.001 Subroutine SGANT calculates this 
quantity by setting up one filamentary dipole on the wire axis and 
another identical dipole on the wire surface. These dipoles lie In 
parallel planes separated by a distance equal to the wire radius. 

In Table 5, dipoles 1 and 2 have terminals at vertices 1 and 2, 
respectively, and they share the middle segment. Again these dipoles 
lie in parallel planes separated by a distance equal to the wire radius. 
For a one-turn planar polygon wire loop, subroutine SGANT would generate 
the data In Table 4 for the diagonal elements Z-}-} and the data in 
Table 5 for the next elements. 
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TABLE 1 

Input and Output Data for Straight Wire 


Input Data 


002 

“ 2.56 

-1.0 

0.0005 



0.001 

1.00 

1.0 

-1.0 0.0 



1 

1 

1 

1 0 

3 

4 

300. 

0 . 

90. 

0. 90. 

45. 

45. 

1 

2 





2 

3 





3 

4 





4 

5 





0 . 

0 . 

-0.250 




0 . 

0 . 

-0.125 

1 

2 3 

4 

0 . 

0 . 

0 . 

• • • 

0 - 

— • # 

0 . 

0 . 

0.125 

1 2 


4 5 

0 . 

0 . 

0.250 


3 


1 . 

1 . 

1 . 




Output Data 





98.18 

0.0095 

82.97 

43.26 



-.091 

0.080 

-0.091 

0.080 

0.224 

-0.096 

0.0 

90.0 

0.0 

1.615 



0.0 

90.0 

0.0 

0.0 

0.0 

0.608 

0.0 

0.0069 

0.0 

0.377 

0.0 

0.370 

45.0 

45.0 

0.0 

0.0 

0.0 

0.239 


TABLE 2 

Input and Output Data for Square Loop 


Input Data 

0.002 2.56 -1.0 0.0005 


0.001 

1.0 

1.0 -1.0 

1 

1 

1 1 

300. 

0.0 

90.0 0.0 

1 

2 


2 

3 


3 

4 


4 

1 

.. 

0.05 

-0.05 

0.0 

0.05 

0.05 

0.0 

-0.05 

0.05 

0.0 

-0.05 

-0.05 

0.0 

1.0 

1.0 

1.0 

Output Data 


73.10 

.243E-4 

62.94 

-.0078 

.0027 

.0057 

0.0 

90.0 

.8066 

0.0 

90.0 

.0002 

.126E-4 

0.0 

.936E-4 

45.0 

45.0 

. 106E-3 


0.0 

0 1 4 

90.0 45. 45. 



1609.8 

.0029 -.0010 -.0056 

.0 

.0 .0 .0 

.0 .810E-4 .0 

.265E-4 .0 .0 


4 
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TABLE 3 

Input and Output Data for Y Antenna 


Input Data 






002 

2.56 

-1.0 

0.0005 



0.001 

1.0 

1.0 

-1.0 

0.0 


1 

1 

1 

1 

0 

2 4 

300. 

0.0 

90.0 

0.0 

90.0 

45. 45. 

1 

2 





2 

3 




,5 

3 

4 




/ 

3 

5 




A 

0.0 

-0.30 

0.0 

1 

(Z \ 2 , 

/ 

0.0 

-0.15 

0.0 

1 

3* 


0.0 

0.0 

0.0 



\ 

0.1 

0.1 

0.0 




-0.1 

0.1 

0.0 




1.0 

1.0 

1.0 




Output Data 





97.88 

0.013 

75.53 

-0.572 



-.124 

0.081 

0.260 

-0.064 

-0.126 

0.070 

0.0 

90.0 

1.535 

0.0 



0.0 

90.0 

0.748 

0.0 

0.0 

0.0 

0.0103 

0.0 

0.487 

0.0 

0.477 

0.0 

45.0 

45.0 

0.360 

0.170 

0.0 

0.0 


TABLE 4 


Self Impedance of Two-Segment V Dipole Shown in Fig. 2 
Radius: a = 0.001 X 


4> 

h = 0.10X 

h = 0.15X 

h = 0.20X 

h = 0.25X 

30° 

60 

90 

120 

150 

180 


1.4 - j 314 
5.3 - j 337 
• 10.4 - j 340 
15.3 - j 338 
18.9 - j 335 
20.1 - j 335 

3.1 - j 186 

11.0 - j 177 

21.1 - j 163 

30.9 - j 151 
37.7 - j 144 

39.9 - j 142 

6.1 - j 61 

21.3 - j 21 

40.0 + j 9 
57.7 + j 28 

69.3 + j 39 

73.1 + j 42 
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Fig. 2. Symmetric center-fed V dipole. 


TABLE 5 


Mutual Impedance Between Overlapping V Dipoles in Fig. 3 

Radi us : a = O.OOl X 


♦ 

h = 0.10X 

msmm 

h = 0.20X 

h = 0.25X 

60° 

90 

120 

150 

180 

-0.96 + j 338 
0.19 + j 322 
3.29 + j 336 
6.61 + j 346 
8.01 + j 349 

-2.08 + j 285 
1.03 + j 276 
8.40 + j 290 
15.61 + j 299 
18.47 + j 301 

-3.45 + j 275 
3.57 + j 271 
17.86 + j 285 
30.00 + j 291 
34.35 + j 292 

- 4.8 + j 298 

10.1 + j 297 
35.3 + j 309 
52.9 + j 309 

58.2 + j 308 



Fig. 3. Overlapping V dipoles share the middle segment. 


Tables 6, 7, and 8 list the mutual impedance for other configura- 
tions. In all these tables, the data apply to two-segment center-fed 
sinusoidal dipoles with identical segment lengths h. 
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TABLE 6 

Mutual Impedance Between Overlapping V Dipoles in Fig. 4 

Radius: a = 0.001 X 


a 

h = 0.10X 

no 


■bb 

30° 

60 

90 

120 

150 

6.74 - j 314 
3.16 - j 291 
0.06 - j 278 
-1.01 - j 256 
-0.48 - j 207 

16.24 - j 167 
7.68 - j 169 
0.31 - j 172 
-2.39 - j 168 
-1.20 - j 146 

32.17 - j 56 
15.47 - j 76 
1.15 - j 92 
-4.47 - j 101 
-2.40 - j 98 

58.7 + j 49.6 

28.8 + j 14.2 
3.5 - j 12.2 

-7.6 - j 35.5 
-4.5 - j 50.7 



Fig. 4. Overlapping V dipoles share the bottom segment in a 
planar Y configuration. 






TABLE 7 

Mutual Impedance Between the Cop! anar- Skew Linear Dipoles in Fig. 5 

Displacement: d = \ 


0 

■BBS 1 

BBSS 

■BBS 

h = 0.25X 

0° 

0.337 + j 1.952 

0.880 + j 4.759 

1.932 + j 9.547 

4.011 + j 17.7 

15 

0.322 + j 1.884 

0.831 + j 4.585 

1.799 + j 9.180 

3.671 + j 17.0 

30 

0.281 + j 1.684 

0.700 + j 4.082 

1.448 + j 8.128 

2.800 + j 15.0 

45 

0.220 + j 1.369 

0.521 + j 3.301 

1.000 + j 6.519 

1.745 + j 11.9 

60 

0.149 + j 0.964 

0.333 + j 2.310 

0.579 + j 4.524 

0.860 + j 8.1 

75 

0.075 + j 0.497 

0.159 + j 1.187 

0.252 + j 2.308 

0.305 + j 4.1 

90 

0.0 + j 0.0 

0.0 + j 0.0 

0.0 + j 0.0 

0.0 + j 0.0 



Fig. 5. Center-fed coplanar-skew linear dipoles. 












TABLE 8 

Mutual Impedance Between the Nonplanar- 

Displacement: 


♦ 

O 

r— 

• 

° 

II 

x: 

masm 

0° 

0.337 + j 1.952 

0.880 + j 4.759 

15 

0.326 + j 1.886 

0.850 + j 4.596 

30 

0.292 + j 1.691 

0.762 + j 4.121 

45 

0.238 + j 1.380 

0.622 + j 3.365 

60 

0.169 + j 0.976 

0.440 + j 2.380 

75 

0.087 + j 0.505 

0.228 + j 1.232 

90 

0.0 + j' 0.0 

0.0 + j 0.0 


Skew Linear Dipoles in Fig. 6 
d = X 

h = 0.20X h = 0.25X 

1.932 + j 9.547 4.011 + j 17.74 

1.867 + j 9.222 3.877 + j 17.14 

1.675 + j 8.269 3.482 + j 15.37 

1.369 + j 6.752 2.850 + j 12.55 

0.969 + j 4.775 2.020 + j 8.88 

0.502+j2.472 1.047 +j 4.60 

0.0 + j 0.0 0.0 + j 0.0 









III. SUMMARY 

This report presents the sinusoidal -reaction computer program for 
thin-wire antennas and scatterers, instructions for the user, typical 
input and output data and mutual -impedance tables for sinusoidal dipoles. 
Appendices list the computer subroutines and explain their functions. 
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APPENDIX 1. Subroutine SORT 

Subroutine SORT, listed in Fig. 7, defines a set of dipole mode 

iMM r f nts ? n the w i re structure - The input data IA, IB, NM, NP, ICJ and 
INM have been defined already. The output data are defined as follows: 

N total number of dipole modes 

11(1) endpoint of dipole I 

12(1) terminal point of dipole I 

13(1) endpoint of dipole I 

JA(l) first segment of dipole I 

JB(l) second segment of dipole I 

MD(J,K) list of dipoles sharing segment J 

ND(J) total number of dipoles sharing segment J 

MAX, MIN extreme values of ND(J) 

At completion of the DO LOOP ending with statement 20, NJK denotes 
the number of segments intersecting at point K, and JSP is a list of 
these segments. In the DO LOOP ending with statement 22, the computer 
sets up the appropriate number MOD of dipoles modes with terminals at 
poink K. 

APPENDIX 2. Subroutine SGANT 

Subroutine SGANT, listed in Fig. 8, calculates the mutual imped- 
ances Z i • and stores them in C(K). The input data for SGANT have been 
defined already. The output data are defined as follows: 

C(K) open-circuit impedance matrix 

CGD(J) cosh yd for segment J 

SGD(j) sinh yd for segment J 

D(J) length of segment J 

ZS surface impedance of the wire 

The surface impedance is calculated just above statement 12. B01 

denotes Jq/Ji where J Q and are the Bessel functions of order zero 
and one with complex argument ZARG. It is assumed that all the wire 
segments have the same radius, conductivity and surface impedance. 

In the DO LOOP ending with statement 20, SGANT calculates the 
segment lengths D(J). DMIN and DMAX denote the lengths of the 
shortest and longest segments. If the wire radius or the segment 
lengths are clearly beyond the range of thin-wire theory, N is set to 
zero at statement 25 followed by RETURN to the main program to abort 
the calculation. 

At statement 30, the program selects a segment K, and a few state- 
ments below this it selects another segment L. K is a segment of test 
dipole I, and L is a segment of expansion mode J. The mutual impedance 
between segments K and L is obtained by calling subroutine GGS or GGMM. 


19 



SUBROUTINE SORT! I A , I B t 1 1 tI2,I3,JA T JB,MDfNDtNM f NP,NtMAX T MlN 

2? ICJtINM) 

DIMENSION JSPI20) 

DIMENSION Il(l)»I2(l)fI3(l)»JA(l)»JB(l) 

DIMENSION I A ( 1) ,IB(1),ND( 1 ) ,MD ( INM ,4 ) 

1=0 

DO 24 K= 1 »N P 
N J K = 0 

DO 20 J = 1 ,NM 
IND=UA(J)-K)*UB(J)-K) 

IF ( IND.NE *0 ) GO TO 20 
N JK=N JK+1 
JSP (NJK )=J 
20 CONTINUE 
MOD =N JK-1 

IF {MOD .LE . 0 ) GO TO 24 
DO 22 IMD= I fMOD 
1 = 1 + 1 

IF ( I *C»T * I C J ) GO TO 22 

I PD = I MD + 1 

JA I = J SP ( IMD ) 

J A { I ) = J A I 
JB I = J SP ( I PD ) 

J B ( I ) = JB I 

I I <i) = ia( jai > 

IF ( I A( JAII .EO.K >1 1 ( I )=IB( JAI ) 

1 2 n ) =k 
I 3 < I ) = I A ( JB I ) 

IF { 1 A( JBI ) .EO*K) I3U ) = IB( JBI ) 

22 CONTINUE 
24 CONTINUE 
N = I 

DO 30 J= 1 >NM 
NO ( J ) =0 
DO 30 K=l»4 
30 MD < J,K )=0 
I I I =N 

IF ( N .GT • I C J ) I 1 1 = I C J 
DO 40 1 = 1, I 1 1 
J = J A ( I ) 

DO 38 L=l,2 
ND { J ) =ND ( j ) + 1 
K = 1 
M = 0 

32 MJK=MD(J,K) 

IF { M JK *NE ,0)00 TO 34 
M = 1 

MD ( J , K ) = I 
34 K = K+ 1 

I F { K • GT .4 ) GO TO 38 
IF {M .EQ • 0 ) GO TO 32 
38 J=JB{ I ) 

40 CONTINUE 
M IN = 1 00 
M AX = 0 

DO 46 J=1,NM 
NO J=ND ( J J 

IF (NDJ .GT .MAX )MAX=NDJ 
4 6 IF (ND J.LT .M IN )M IN =ND J 
RETURN 
END 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0018' 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 


Fig. 7. Subroutine SORT 
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SUBROUTINE SGANT (I A , IB , INM , 1NT , I SC , 1 1 , I 2 , I 3 * JA, J B ,MD , N , ND , NM , NP 
2, AM ,BM » C iCGO ,CMM ,D ,E P2 , E P3 , E T A ,F HZ ,G AM , SGD,X , Y , Z ,Z L D, Z S ) 

COMPLEX ZG,ZH,ZS,EGD ,GD ,CGDS , SGO S , SGDT, BO 1 

C0MPLE X Pllt P12 f P21 ,P22,Q11 ,012 ,021 , 0 22, EP2, EP , ETA, GAM, EP3 

COMPLEX EPS ILA,CWEA, BETA t ZARG 

COMPLEX P(2,2)»Q(2,2),CGD(1 >,SG0(1),C( 1 J ,.ZLD( 1 ) 

D I MEM SION XC 1 J ,Y( 1 J ,Z( 1) ,0 ( 1 ) ,1 A( 1 ) ,1 B( 1 ) ,MD( INM, 4) 

D IMENS ION 11(1), 12(1), 13(1), UA(1),JB( 1),ND( 1),ISC(1) 

DATA EO , TP ,U0/8.854E -12 ,6.283 18 , 1 .25666-6/ 

2 FURMATC 3X , 'AM = * ,E 1 0 . 3 , 3X , 'DM AX = * ,E 10. 3, 3X , f DM I N = %E10.3> 

E P=E P3 

ICC = (N*N+N )/2 
DO 10 1=1, ICC 
]0 c<n = (,o,.o) 

ZS= ( .0, .0) 

IF (CMM.LE • 0 . ) GO TO 12 
OMEGA = TP*FHZ 

EPS I LA = CM PLX(E0,-CMM*1.E 6/OMEGA) 

CWEA=(*0,l.)*OMEGA*EPSILA 

BE TA=OMEG ASSORT (UO )*CSQRT (E PSILA-EP) 

Z ARG= BE T A=*AM 

CALL CBE S ( ZARG , BO 1 ) 

ZS=BETA*B01/CWE A 
]2 ZH=ZS/( TP*AM*GAM ) 

DM IN = 1 .E 30 
DM AX= .0 
DO 20 J= 1 ,NM 
K= I A ( J ) 

D ( j ) = SORT ( (XCK )— X « L) ) **2 + ( Y (K ) -Y ( L ) ) ** 2+( Z ( K> -Z ( L ) ) **2> 

IF (0( J) .LT.DMIN )DMIN=D( J) 

IF (D( J) .GT.DMAX )DMAX=D( J) 

EGD=CEXP(GAM*D( J) ) 

CGD(J) = (EGD + l,/EGD)/2* 

20 SGDU)=(EGD-1./EGD )/2. 

IF (DMIN,LT.2.*AM )G0 TO 25 
IF (CABS(GAM*AM ) .GT.0.06IG0 TO 25 
IF ( CABS ( GAM *DM AX).GT.3« )G0 TO 25 
IF ( AM .GT .0. )G0 TO 30 

25 * = 0 

WRITE (6, 2) AM , DM AX, DM IN 
RETURN 

30 DO 200 K = 1 ,NM 
NDK=ND( K ) 

K A= I A ( K ) 

KB= I B ( K ) 

DK=D ( K) 

CGDS=CGD(K) 

SGD S= SGO ( K ) 

DO 200 L=1,NM 
ND L = ND ( L ) 

LA = I A ( L ) 

LB= I B ( L ) 

DL = D( L) 

SGD T = SGD ( L ) 

N I L = 0 

DO 200 1 1 = 1 , NDK 
I =MD ( K , 1 1 ) 

M M = ( I - 1 ) *N - ( I * I — I ) / 2 
F 1 = 1 . 

IF (KB.EQ.I2 ( I ) ) GO TO 36 
IF ( KB.EQ , 1 1 ( I ) )F 1 =“!• 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 
00 13 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
00 21 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 
004 8 

0049 

0050 
00 51 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 


Fig. 8a. Subroutine SGANT 
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I S= 1 0063' 

GO TO 40 0064 

36 IF(KA*EQ.I3(I MFI--1* 0065 

15=2 0066 

40 DO 200 JJ=1 t NDL 0067 

J=MD(LtJJ) 0068 

MMM=MM+J 0069 

IF ( I *GT * J ) GO TO 200 0070 

F J= 1 • 0071 

IF ( LB *E Q * 12 ( J ) )G0 TO 46 0072 

IF ( LB.EQ.I 1 ( J ) )F J=-l • 0073 

JS=1 0074 

GO TO 50 0075 

46 IF ( LA*EQ.I3(J) JF J=-l • 0076 

JS=2 0077 

50 IF (NIL *NE .0 ) GO TO 168 0078 

N I L= 1 0079 

IF (K.FQ.L ) GO TO 120 0080 

IND = ( LA-KA )*( LB-KA ) * ( L A-KB )*(LB-KB> 0081 

tel IND.EO.OJGO TO 80 0082 

C SEGMENTS K AND L SHARE NO POINTS 0083 

CALL GGS ( X ( K A ) , V { KA ) , Z ( K A ) ? X ( KB ) , Y ( K B ) , Z { KB J ? X ( LA ) , Y( L A ) ,Z ( LA ) 00 84 

2»X(LB),YILB),Z( LB) t AM f DK,CGDS,SGDS,DL,SGDT r INT,ETA f GAM 0085 

3« P( 1 »1 ) »P(1 t2) «P(2 yl ) tP(2 f 2 ) ) 00 86 

GO TO 168 0087 

C SEGMENTS K AND L, SHARE ONE POINT (THEY INTERSECT) 0088 

80 KG=0 ' 0089 

JM=KB 0090 

JC=KA 0091 

KF=1 0092 

INO=(KB-LA )*( KB-LB) 0093 

IF ( I ND *NE *0 } GO TO 82 0094 

JC= KB 0095 

KF=-1 0096 

JM = K A 0097 

K G= 3 0098 

82 LG= 3 0099 

J P= LA 0100 

LF=-1 0101 

IF (LB.EQ.JC)GO TO 83 0102 

JP=LB 0103 

LF = 1 0104 

LG=0 0105 

83 SGN = KF*LF 0106 

CPSI= U X( JP )-X( JC) >*(X( JM)-X { JC )) +( YCJPJ-Yl JC) )*( Y( JM) -Y< JC) ) 0107 

2+( U JP )-Z ( JC ) )*( Z ( JM }-Z ( JC) ) )/{DK*DL ) 010 8 

CALL GGMM ( .0 ,DK f . 0 t D L t AM T C GD S , SGO S , SGD T f CP SI , ETA, GAM 0109 

2» 0( 1 y 1 ) ,Q( 1 v ? ) v 0(2 »l)’yQ(2f2 )) OHO 

DO 98 KK = 1 ,2 0111 

KP = I ARS ( KK-KG ) 0112 

DU 9 8 L L = 1 » 2 0113 

L P = I A B S ( LL-LG) 0114 

P(KP,LP> =SGN*Q(KK t LL) 0115 

98 CONTINUE 0116 

GO TU 168 0117 

C K=L (SELF REACTION OF SEGMENT K) 0118 

120 Q11=(.0».0> 0119 

012=1 • 0 y • 0 ) 0120 

IF ( CMM * LE . 0 * ) GO TO 150 0121 

GD = GAM*DK 0122 

ZG = ZH/( SGD S **2 ) 0123 

Q11=ZG*(SGDS*CGDS-GD ) /2 . 0124 

Fig. 8b. Subroutine SGANT 
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012«ZG*(GD*CGDS-SGDS>/2. 0125 

150 I SCK=I SCI K ) 0126 

PllM.0,.0) 0127 

P12=(.0,.0) 0128 

IF ( I SCK.EO .0 ) GO TO 155 0129 

IF (BM.LE .AMIGO TO 155 0130 

CALL DSHBLLIAM,BMfDK,CG0S»SG0SfEP2»EP,ETAtGAM,Pll,P12) 0131 

155 Q11=P11+Q11 0132 

Q12=P12+Q12 0133 

CALL GGMM(.0,DK,.0,DK,AM,CGDS,SGDS,SGDS,1. 0134 

2 t E T A * GAM »P11»P12*P21*P22) 0135 

011=P11+011 0136 

012=P12+01 2 0137 

P(1»1)=Q11 0138 

P( 1 *2 ) = Q 1 2 0139 

P ( 2 * 1 J =Q12 0140 

P ( 2 » 2 ) = Q1 1 0141 

IF IKA.NE .LAJGO TO 160 0142 

GO TO 168 0143 

160 P ( 1 * 1 ) = — Q 1 2 0144 

P ( 1 t 2 ) =-Ql 1 0145 

P ( 2 » 1 ) =-Q 11 0146 

P(2,2)=-012 0147 

168 C(MMM)=C(MMM )+FI*FJ*P(IS, JS ) 0148 

200 CONTINUE 0149 

DO 220 I = 1 *N 0150 

I J»Cl-U*N-U*I-n/2 + I 0151 

J 1 = J A < I ) 0152 

IF ( 12 ( I ) .EQ.IBI J1 )> J1 = J1+NM 0153 

J2=J8(I) 0154 

IF (12(1) .EQ.IB( J2) JJ2=J2+NM 0155 

220 CIIJ)=C(IJ)+ZLD(J1)+ZLD(J2) 0156 

RETURN 0157 

ENO 0158 


Fig. 8c. Subroutine SGANT 
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In statement 168, this impedance is lumped into C(MMM). The mutual 
Impedance Z. . between dipoles I and J is the sum of four segment- 
segment impedances. 

In SGANT, segment K has endpoints KA and KB, and segment L has 

endpoints LA and LB. It is convenient to think of KA and KB as points 

1 and 2 on segment K, and LA and LB as points 1 and 2 on L. Now we 
define four segment-segment impedances P(IS,JS). The first subscript 
IS refers to the terminal point on segment K, and the second subscript 

JS refers to the terminal point on L. Thus IS = 1 or 2 if dipole I 

has its terminal point 12(1) at KA (point 1) or KB (point 2), respec- 
tively. Similarly, JS = 1 or 2 if mode J has its terminal point I2(J) 
at LA or LB. The impedances P(IS,JS) are defined with the following 
reference directions for current flow: from point 1 toward point 2 on 

each segment. If dipole I has this same reference direction on 
segment K, we set FI = 1; otherwise FI = -1. Similarly FJ = 1 or -1 
in accordance with the reference direction for mode J on segment L. 

In statement 168, P(IS,JS) is multiplied by FI and FJ before its 
contribution is added to Z^j. 

Subroutine GGMM calculates the impedances Q(KK,LL) which are like 
the P(IS,JS) but have different conventions for reference directions 
and subscript meaning. The transformation from the Q impedances to the 
P impedances Is accomplished in the DO LOOP ending with statement 98. 

If the wire has finite conductivity, the appropriate modifi- 
cation is applied to the impedance matrix just above statement 150. 

(See Eqs . 27 through 29 in Reference 1.) The terms arising from the 
dielectric shell on an insulated segment are obtained from subroutine 
DSHELL just above statement 155. Finally, the lumped loads ZLD are 
added to the diagonal elements of the impedance matrix In statement 220. 

The inpedance matrix could be calculated in a different order as 
follows. Select modes I and J, calculate ZIJ, and then increment I 
or J. Instead, SGANT selects segments K and L, calculates ZKL, adds 
ZKL to all the appropriate elements ZIJ, and then Increments K or L. 

This minimizes the calls to GGS and GGMM and presumably improves the 
computational efficiency. 

K is a segment of test dipole I, and L is a segment of expansion 
mode J. When the segment numbers K and L are equal, SGANT calls GGMM 
to obtain the mutual impedance between two filamentary electric mono- 
poles. These monopoles are parallel and have the same length. Mono- 
pole K is positioned on the axis of the wire segment, and monopole L 
Is on the surface of the same wire segment. Thus, the displacement is 
equal to the wire radius. The two monopoles are side-by-side with no 
stagger. 

When segments K and L intersect, SGANT again calls GGMM for the 
mutual Impedance between the two filamentary monopoles. Monopole K is 
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situated on the axis of wire segment K, and monopole L is on the surface 
of wire segment L. The axes of segments K and L define a plane P, and 
monopole K lies in this plane. Monopole L is parallel with plane P and 
is displaced from it by a distance equal to the wire radius. 

APPENDIX 3. Subroutine CBES 

Subroutine CBES, listed in Fig. 9, calculates the quantity B01 = 
Jq(z)/Ji(z) where z is complex and J Q and denote the Bessel func- 
tions of order zero and one. 

APPENDIX 4. Subroutine DSHELL 

Subroutine DSHELL, listed in Fig. 10, calculates the mutual im- 
pedance term contributed by the dielectric insulation on the surface 
of a thin wire. This subroutine uses Eq. 35 of Reference 1. 

APPENDIX 5. Subroutine GGS 

Subroutine GGS, listed in Fig. 11, calculates the mutual impedance 
between two filamentary monopoles with sinusoidal current distributions. 
(The dipole-dipole mutual impedance in Eq. 20 of Reference 1 is the sum 
of four monopole -monopole mutual impedances.) The endpoints of the 
axial test monopole s are (XA,YA,ZA) and (XB,YB,ZB), and the endpoints 
of the expansion monopole t are (X1,Y1,Z1) and (X2,Y2,Z2). DS and DT 
denote the lengths of monopoles s and t, respectively. CAS, CBS and 
CGS are the direction cosines of monopole s, and CA, CB and CG are the 
direction cosines of monopole t. 

If INT = 0, GGS calls GGMM for the closed- form impedance calcula- 
tions. Otherwise GGS calculates the mutual Impedance via Simpson's- 
rule integration with the following nunfoer of sample points: IP s 

INT + 1. If the monopoles are parallel with small displacement, GGS 
calls GGMM to avoid the difficulties of numerical Integration. 

For the fields of the test monopole, GGS uses Eqs. 75 and 76 of 
Reference 1. The current distribution on the expansion monopole is 
given by Eq. 74 of Reference 1. With an origin at (X1,Y1,Z1), the 
coordinate T measures distance along the expansion monopole. Thus T 
is the integration variable. 

Let the coordinate s measure distance along the test monopole with 
origin at (XA,YA,ZA). From any point T on monopole t, construct a 
line to the test monopole such that the line is perpendicular to the 
test monopole. SZ denotes the s coordinate of the intersection of this 
line with the test monopole. The length of the line is the radial 
coordinate p, and RS denotes p 2 . R1 and R2 are the distances from 
(XA,YA,ZA) and (XB,YB,ZB) to the point T. Cl is the current at T for 
the mode with terminals at (X1,Y1,Z1), and C2 is the current at T for 
the other mode with terminals at (X2,Y2,Z2). C denotes the Simpson's- 
rule weighting coefficient. 
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SUBROUTINE CBES(Z,B01) OOOl 

COM RLE X ARG,CC,CS,EX 0002 

COMPLEX BOI ,Z,TERMJ,TERMN ,MZ24, JN { 2) 0003 

DATA PI /3* 14 159/ 0004 

IF ( CABS < Z J.GE.12.0) GO TO 10 0005 

FACTORS. 0 0006 

TE RMN = ( 0. ,0. J 0007 

MZ24=-0.25*Z*Z 0008 

TERMJ=( 1. 0,0.0) 0009 

DO 1 NP= 1 , 2 0010 

N=NP-1 0011 

JN(NP)=TERMJ 0012 

M = 0 0013 

2 M=M+ 1 0014 

TERMJ = TERMJ*MZ24/FL0AT(M*(N+M ) ) 0015 

JN (N P ) = JN IN P ) +TERMJ 0016 

IF (NP.NE • 1 ) GO TO 3 0017 

FACTO R=F ACTOR + 1 •O/F LOAT(M) 0018 

TE RMN = TE RMN+TE RMJ*F ACTOR 0019 

3 ERROR = CAB$(TERMJ) 0020 

IF (ERROR.GT *1 . 0E-10 ) GO TO 2 0021 

1 TERM J = 0 .5 *Z 0022 

BO 1 = JN ( I )/JN(2 J 0023 

RETURN 0024 

10 Y = A IM AG I Z ) 0025 

IF {ABS(Y) .GT *20. )G0 TO 20 0026 

ARG= ( *0,1. )*Z 0027 

E X = CE XP [ ARG ) 0028 

CC=EX+1./EX 0029 

CS = ( .0,-1 • >* (EX-1. /EX) 0030 

B01 = < CS+CC)/(CS-CC) 0031 

RETURN 0032 

20 801=1*0,-1.) 0033 

IF (Y.LT.O. >BOl=( *0,1.) 0034 

RETURN 0035 

END 0036 


Fig, 9. Subroutine CBES 
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SUBROUTINE D SHE LL ( AM , BM ,DK , CGD S, $GO S »E P2, EP , E TA , GAM , P 1 1 , PI 2 ) 
COMPLEX CGDS,SGDS,EP2,EP,ETA,GAM ,P 1 1,P12 , GO, CS T 
DATA PI/3. 14159/ 

GD = GAM*DK 

CST=(E P2-EP )*ETA*ALOG( BM/AM ) /(4 .*PI *E P2*SGDS*SGDS) 
Pll=-CST*(GD+$GOS*CGD$) 

P12=CST*(GO ♦CGD S+SGDS ) 

RETURN 

END 
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Fig. 10. Subroutine DSHELL. 
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SUBROUTINE GGS ( XA , YA ,ZA t XB ,YB»ZB,X 1, V 1 ,Z 1 »X 2* Y2 ,Z2 t AM 
2?DS»CGDSt$GDS ,DT,SGDT t INT t ETA,GAM t Pll f P12,P21fP22) 

COMPLEX Pll ,P12 ,P21,P22 ,E JA t E JB,E Jl»E J2, ETA, GAM, Cl » C2tCST 
COMPLEX EGO »CGO S » SGD $ ? SCOT ,ERl»ER2 f ETl,E T2 
DATA F P / 12 .5 6637/ 

CA=(X2-X1)/DT 
CB=< Y2-Y1 )/DT 
CG=(Z2-Z1)/DT 
CAS=(XB-XA)/DS 
CBS=IYB~YA)/DS 
CGS=(ZB-ZA)/DS 
CC=CA*CAS+CB*CBS+CG*CGS 
IE ( ABS( CC ) .GT • *997 )GO TO 200 
20 SZ= ( Xl-XA ) *CAS+(Yl-YA)*CBS + ( Z 1-Z A) *CGS 
IF ( INT .LE .0)G0 TO 300 
IN S = 2 * I INT/2 ) 

IF (INS.LT.2)INS*2 
I P= IN S+ 1 
DE LT=DT /INS 
T = *0 

DSZ=CC*DELT 
Pll=< .0,.0) 

P12=( .0».0) 

P2l=<.0,.0) 

P22 = I .0,.0) 

AMS = AM *AM 
SGN =-l • 

00 100 I N= 1 1 1 P 

ZZ1=SZ 

ZZ2=SZ-DS 

XXZ=X1+T*CA-XA-SZ*CAS 

YYZ-Y1+T*CB-YA-SZ*CBS 

ZZZ=Z1+T*CG-ZA-SZ*CG$ 

RS=XX2**2+YYZ**2+ZZZ**2 
R1=SQRT<R$+ZZ1**2 ) 

E JA=CEXP(-GAM*Rl) 

E J1=E JA/R1 
R2=$QRT (RS+ZZ2**2 ) 

EJB=CEXP(-GAM*R2 ) 

E J2 = E JB/R2 

ER1=E JA*SGDS+ZZ 1 *E J 1 *CGD S-Z Z2 *E J2 
ER2=-EJB*SGDS+ZZ2*E J2*CGDS-ZZ 1*E J1 
F AC= *0 

IF !RS.GT.AMS)FAC*<CA*XXZ+CB*YYZ+CG*ZZZ)/RS 
ET1 = CC*IE J2-EJl*CGDS)+FAC*ERl 
ET2=CC*(EJ 1-E J2*CGDS )+FAC*ER2 
C = 3 • + SGN 

IF ( IN .E Q • l .OR. IN .E Q» IP )C c l • 

EGD = CEXPf GAM*(DT-T }) 

C1=C*(EG0-1 ./EGOJ/2. 

E GD=CE XP { GAM*T ) 

C2 = C*(EG0-l./EGD)/2* 

PI 1 = P1 1+E Tl*Ci 
P12=P12+ET1*C2 
P21=P21+ET2*C1 
P22 = P22+ET2 *C2 
T= T +0E LT 
SZ=$Z+DSZ 
100 SGN =- SGN 

CST=-ETA*0ELT/(3**F P*SGOS*SGDT ) 

P11=CST*P11 
P12=CST*P1 2 
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Fig. 11a, Subroutine GGS 
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P2 1*CST *P2 1 
P22=CST*P22 
RETURN 

200 SZ1=(X1-XA)*CAS+(Y1-YA)*CBS+(Z1-ZA)*CGS 

RH 1 = SQRT ( (X1-XA-SZ1*CAS )**2+{ Y 1-YA-SZ l*CB$)**2+( Z l-ZA-SZ 1*CGS)**2> 
S Z 2= SZ 1 +DT *CC 

RH2 - SO RT ( ( X2-XA-SZ2*CAS )**2 +( Y2-YA-SZ 2*CB5)**2+( Z2-ZA-SZ 2*CGS)**2) 
DDD= (RH1+RH2 )/2* 

IF (ODD .GT .20**AM .AND. INT.GT.O)GO TO 20 
IF (DDD.LT.AM )DDD=AM 

CALL GGMM ( .0 ,D S , S Z 1 t SZ2 ,DDD ,CGDS,SGD S, SGDT, 1. 

2,ETA,GAM,P11,P12,P21 ,P22) 

RETURN 

300 $S*SQRT( l.-CC*CC) 

CAD=(CGS*CB-CBS*CG)/SS 

CBD = ( CAS*CG-CGS*CA ) / SS 

CGD=(CB5*CA-CAS*CB)/SS 

DK= ( Xi-XA) *CAD+( Y1-YA)*CBD+(Z WA)*CGD 

DK=AB$ ( DK ) 

IF (DK.LT .AM )D K=AM 

XZ=XA+SZ*CAS 

Y Z =Y A + SZ *CB S 

ZZ=Z A+SZ *CGS 

X P 1 = X 1-D K*CAD 

YP1=Y1-DK*CBD 

ZP1 = Z1-DK*CGD 

CAP=CBS*CGD-CGS*CBD 

CBP=CGS*CAD-CAS*CGD 

CGP=CAS*CBD-CBS*CAD 

P1=CAP*( XP1-XZ)+CBP*(YP1-YZ ) +CGP*t ZPWZ ) 

T1=P1/ SS 
S1=T1*CC-SZ 

CALL GGMM (SI ,S1+DS ,T1 ,T1+DT ,DK,CGO S, SGDS, SGDT, CC, ET A, GAM 
2tPlltP12,P21,P22) 

RETURN 

END 
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Fig. lib. Subroutine 6GS 
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Below statement 300, GGS performs some analytic geometry In prepara- 
tion for calling GGMM. The remaining part of this Appendix concerns 
this last part of subroutine GGS. 

Let s denote a unit vector in the direction from (XA,YA,ZA) toward 
(XB,YB,ZB|. Also let t denote a unit vector from (X1,Y1,Z1) toward 
(X2,Y2,Z2). Then s • t = cos 0 = CC where e is the angle formed by the 
axes of the two monopoles. Let monopole s lie in one plane P and 
monopole t lie in another parallel plane P t . ^CAD, CBD and CGD are the 
direction cosines of the unit vector d = t x s / sin 0 which is perpen- 
dicular to both planes. To obtain the distance DK between the two 
planes, we construct a vector R_n from (XA,YA,ZA) to (X1,Y1,Z1) and 
take DK = R n • d. 

Construct a line from (X1,Y1,Z1) to the test monopole, such that 
the line is perpendicular to the test monopole. SZ denotes the s co- 
ordinate of the intersection of this line with the test monopole, and the 
cartesian coordinates of this intersection are XZ, YZ and ZZ. The 
direction cosines of s x £ are CAP, CBP and CGP. 

From the point (X1,Y1,Z1) in plane P-£, construct a perpendicular 
line to the point (XP1,YP1,ZP1) in the plane P s . This line is parallel 
with a and has length DK. Let II represent a vector from (XZ,YZ,ZZ) 
to (XP1,YP1,ZP1). PI denotes R • (?xJ). SI and T1 are defined in 
the next Appendix. 

APPENDIX 6. Subroutine GGMM 

Subroutine GGMM calculates the mutual impedance between two fila- 
mentary monopoles with sinusoidal current distributions. The dipole- 
dipole mutual impedance in Eq. 20 of Reference 1 is the sum of four 
monopole-monopole mutual impedances. The monopole impedances are cal- 
culated by GGS with Simpson's rule or by GGMM with closed-form expres- 
sions in terms of exponential Integrals. 

To explain the input data for GGMM, reference is made to Fig. 12. 
Subroutine GGMM is listed in Fig. 13. If the monopoles are parallel, 
let the z axis be parallel with both monopoles. The coordinate origin 
may be selected arbitrarily. SI and S2 denote the z coordinates of the 
endpoints of the test monopole, T1 and T2 are the z coordinates of the 
endpoints of the expansion monopole, and D is the perpendicular distance 
(displacement) between the monopoles. The mutual impedance of parallel 
monopoles Is calculated in the last part of GGMM below statement 110. 

For skew monopoles, let the test monopole s lie in the xy plane 
and the expansion monopole t in the plane z = D. (D is the perpendi- 
cular distance between the parallel planes.) If the monopoles are 
viewed along a line of sight parallel with the z axis as in Fig. 12, 
the extended axes of the two monopoles will appear to intersect at a 
point on the xy plane. Let s measure the distance along the axis of 
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L< — D — >1 




SKEW MONOPOLES 


Fig. 12. Coordinates for parallel and skew monopoles 
in subroutine GGMM. 


31 



SUBRUUT INE GGMM ( S 1 » S2 , T 1 , T2 ,D , CGO S , S GDI , SGD2 , CP SI , E TA,GAM 


OOUBLP PRECI SION Rl»R2fDPO»SIS»TSl*TS2fSTlfST2»COtBOfCPSSfSK 
2»TL1,TL2,TD1,TD2 > SD I ,D PS I ,DD , ZD 
COMPLEX CGDS,SGDS,SGDT,SGD1 ,SGD2,E TA,GAM f P 11 , P 12, P2 1, P22 
COMPLEX CST »EB,EC,EK f ELrEKL,EGZI ,E S1,ES2 , ET1 ,ET2, EXP A, EXPB 
COMPLEX E ( 2 ,2 ) ,F ( 2 ,2 J 


COMPLEX EGZ(2,2),GM(2),GP(2) 
DATA PI/3. 14159/ 


OSG=D*D 


S GO S = S GO 1 

IE (S2.LT.S1 ) SGO $=-SGD 1 
SGO T = SG02 

IE (T2.LT.T1 5 SCOT =-SGD2 

IF (ABStCPSI ).GT..997)GO TO 110 

ES1=CEXP(GAM*S1) 

ES2 = CEXP(GAM*S2) 
ET1=CEXP(GAM*T1 ) 

E T2 = CE XP ( GAM *T 2 ) 


DD=D 

OPSI=CPSI 
TO 1 = T 1 
TD2 = T2 

CPSS=DPSI*DPSI 
CD=00/D SO RT ( l.DO-CPSS) 

C=CD 

BD = CD*DPSI 
B = BD 

EB = CE XP { GAM *CM PLX ( .0,8) ) 

EC=CEXP(GAM*CMPLX( .0,0 ) 

00 10 K= 1 * 2 
DO 10 L= 1 , 2 
10 E<K,L)=(.0,.0) 

TS 1=TD 1 *T0 1 
T S2 = Tf)2*TD2 
DPQ = DD*DD 
S I = S 1 

00 100 1=1,2 

F I = < -1 )**I 

SO I = S I 

S I S = SD I *S0 I 

ST1=2.*SDI*TD1*0PSI 

ST2=2 • *SD I *TD2 *0 PS I 

R 1=0 SORT (DPQ+SIS+TS1-ST1 ) 

R2=DS0RT ( DPG+S I S+T52-ST2 ) 

EK = E B 

DO 50 K=l,2 

FK=(-1) 

SK=E K*SD I 
E L=EC 

DO 40 L=1 ,2 
FL=(-1>**L 
EK L=EK *E L 
XX=FK*BD+FL*CD 
TL1=E L*TD 1 
TL2=F L*TD2 
RR 1 =R1 + SK + TL1 
RR2=R2+SK+TL2 

CALL E XPj <GAM*CMPLX{RR1 ,-XX ) ,G AM *CMP LX ( RR2, -XX ) , EX P A) 
CALL EXP J ( GAM ♦CMP LX { RR 1 ?XX ) ,G AM *CM P LX ( RR2 ,XX } , EXPB) 
E(K,L)»E <K,L)+FI*(EXPA*EKL+EXPB7EKL> 

40 ELM. /EC 
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Fig. 13a. Subroutine GGMM 
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50 EK*1./EB 0063 

ZD = SDI*0PSI 0064 

ZC=ZD 0065 

EGZI=CEXP(GAM*ZC) 0066 

RR 1=R1 +Z0-TD 1 0067 

RR2=R2+ZD-TD2 0068 

CALL EXPJ(GAM*RR1 »GAM *RR2 »E XP B ) 0069 

RR1=R1-ZD+TD1 007 0 

RR2 =R2 -ZD +TD2 ° 071 

CALL EXPJ(GAM*RR1 ,GAM*RR2 ,EXPA) 0072 

F ( I ,1)=2.*SG0S*EXPA/EGZI 0073 

F ( I ,2 )=2.*SGDS*EXPB*EGZI 0074 

100 S I = S2 0075 

CST=ETA/(16.*PI*SGDS*SGDT) 0076 

P 1 1 = CST* < ( F (1 ,1 )+E (2,2 )*ES2-E ( 1,2)/ES2)*ET2 0077 

A +(-F ( 1,2) — E (2, 1 ) *E S2+E ( 1,1) /ES2)/ET2) 0078 

P12=CST*( (-F (1 ,1 )-E (2 ,2 )*ES2+E ( 1 , 2 ) /E S2 ) *ET1 0079 

B +( F(1,2)+E(2,1)*ES2-E(1,1)/ES2)/ET1) 0080 

P21 = CSTM (-F !2 ,1 )-6 (2,2 )*ES1+E ( 1 , 2 ) /E SI ) *ET2 0081 

C +( F(2,2)+E(2,1)*ES1-E(1,1)/ES1)/ET2> 0082 

P22 = CST* ( ( F (2 tl )+E (2,2 )*ES 1-E ( 1,21/E SI ) *ET1 0083 

D +<-F(2,2)-E(2,l)*ESl+E( 1 , 1 ) /E SI ) /E Tl ) 0084 

RETURN 0085 

110 IF (CPSI .LT.O. )G0 TO 120 0086 

T A = T 1 0087 

TB = T2 0088 

GO TO 130 0089 

120 T A = -T 1 0090 

TB=-T2 0091 

SGDT = -SGDT 0092 

130 SI=S1 0093 

DO 150 1=1,2 0094 

T J = TA 0095 

DO 140 J=1 ,2 0096 

Z I J = T J - S I 0097 

R=SORT (DSQ+ZIJ*ZIJ) 0098 

W - R + Z I J 0099 

IF (ZIJ.LT .0.)W=DSQ/(R-ZIJ) 0100 

V =R-Z I J ° 101 

IF (ZlJ.GT.O.)V=DSQ/(R+ZIJ) 0102 

IF ( J.EQ.l )V1=V 0103 

IF ( J.EQ.1)W1=W 0104 

EGZ ( I , J ) =CEXP(GAM*ZIJ ) 0105 

140 T J = TB 0106 

CALL E XPJ (GAM*V1 ,GAM*V ,GP (I )) 0107 

CALL EXPJ (GAM*W 1 ,GAM ,GM ( I ) ) 0108 

150 S I = S2 0109 

CST=-ETA/(8.*PI*SGDS*SGDT) 0110 

PI 1 = CST*< GM (2 )*EGZ (2 ,2 >+GP(2 ) /EGZ (2, 2) 0111 

2-CGDS*(GM(l)*EGZ(l,2)+GP(l)/EGZ(l,2))) 0112 

P12=CST*(-GM (2 )*EGZ <2 , 1 )-GP (2 ) /EGZ(2, 1) 0113 

2 + CGD S* ( GM ( 1 )*EGZ(1,1 J+GPI 1 )/EGZ( 1,1)) ) 0114 

P21=CST*(GM ( 1 ) *EGZ ( 1 ,2 )+GP( 1 ) /EGZ (1,2) 0115 

2-CGDS*(GM<2 )*EGZ(2,2 )+GP(2 )/EGZ(2,2)) ) 0116 

P22=CST*(-GM(1 )#E GZ ( 1 » 1 ) — GP ( 1 ) /EGZ ( 1, 1 ) 0117 

2+CGDS*( GM (2 )*EGZ ( 2, 1 )+GP(2 )/EGZ (2, 1 ) ) ) 0118 

RETURN 0119 

END 0120 


Fig. 13b. Subroutine GGMM 
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the test monopole with origin at the apparent intersection. SI and S2 
denote the s coordinates of the endpoints of the test monopole. Simi- 
larly, let t measure distance along the axis of the expansion monopole 
with origin at the apparent intersection. T1 and T2 denote the t co- 
ordinates of the endpoints of the expansion monopole. Let s and £ be 
unit vectors^ parallel with the positive s and t axes, respectively. 

Then CPSI = s • t = cos The monopole lengths are d s and dt, and 
the remaining input data are defined as follows: 

CGDS cosh yd s 

SGD1 sinh yds 

SGD2 sinh ydt 

GGMM calls EXPJ for the exponential integrals. 

The output data from GGMM are the impedances Pll, P12, P21, and 
P22. In defining these impedances, the reference direction Is from 
SI to S2 for the current on monopole s, and from T1 to T2 for the 
current on monopole t. In the impedance P-,- -j , the first subscript is 
1 or 2 if the test dipole has terminals at si or S2 on monopole s. 

The second subscript is 1 or 2 if the expansion dipole has terminals 
at T1 or T2 on monopole t. The endpoint coordinates SI, S2, T1 and T2 
may be positive or negative. The monopole lengths dc and d t are 
assumed positive in defining the input data CGDS, SGD1 and SGD2. 

For parallel monopoles, CPSI = 1 or -1. SI, S2, T1 and T2 are 
cartesian coordinates for parallel monopoles and spherical coordinates 
for skew monopoles. For skew monopoles, the radial coordinates SI, S2, 

T1 and T2 tend to infinity as the angle if> tends to zero or ir. Therefore, 
if the monopoles are within 4.5° of being parallel, they are approxi- 
mated by parallel dipoles. 

APPENDIX 7. Subroutine EXPJ 

Subroutine EXPJ, listed in Fig. 14, evaluates the exponential 
integral defined as follows: 

f V2 -v . 

(2) W12 ■ I » E,(V1) - E,(V2) + j 2n, 

Jvi V 


where the integration path is the straight line from VI to V2 on the 
complex v plane and 

(3) E,(z) = f Sllit 
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SUBROUTINE E X PJ ( V 1 , V2 ,W 12 ) 

COMPLEX EC,E15 ,S,T,UC,VC,V1 ,V2,W12,Z 
0 IMEN SION V (21 ) ,W (21 ) ,D< 16),E ( 16) 

DATA V/ 0*222 84 667E 00, 

20 • 1 1 889 32 IE 0 1 ,0 .2992 7 363E 01 , 0 .5 775 14 36 E 01 ,0.98 374674E 01, 
20.15982874E 02 , 0 .9 33078 1 2E -0 1 , 0.49 269 174£ 00, 0. 1 21 5 59 54E 01, 

2 0.22 699495E 0 1 ,0 . 366 7622 7E 0 1 , 0 .54 253366 E 01 ,0 .75659 1 62 E 01, 
20.1012022RE 02 ,0 . 1 31 302 82E 02 , 0. 166544 08E 02, 0. 2077 64 79 E 02, 
20.25 623894 E 02,0.314075196 02 , 0 . 385 3068 3E 02 , 0.48026086E 02/ 
DATA W/ 0.45896460E 00, 

20 .4 1 700083E 00 ,0 . 1 1 337 338E 00, 0. 10399 197 E-01 ,0. 26101 720E-03, 

20 . 89854 79 IE-06 ,0 .2 1 82 34 87E 00,0.34 2210178 00, 0. 26302758 E 00, 

20 • 1 2 64 25 82 E 00 ,0 .402 06 865E -0 1 , 0 . 85 6 38778 E-02 , 0. 121 24 3 61 E-02, 

20.1 116744 OE- 03 ,0.64 5 992 67E -05 ,0. 2226 3169 8-06,0.4227 4304E-0 8, 
20.39 2 18973E-10, 0.145 65 152E-12,0.14830270E-15,0. 160059 49 E- 19 / 

DATA 0 / 0.22495842E 02, 

2 0.74411568E 02,-0.414315766 03,-0.787543396 02, 0.11254744E 02, 
2 0.1 602 1 7 61E 03 ,-0.23862 195E 03 , -0.5 0094687E 03 , -0. 6848 7 854E 02, 
2 0.1 225477 8E 02 , -0 . 1 0 1 61976E 02 ,-0.472195918 01, 0.79 729681E 01, 
2-0.210695 74E 02, 0.22046490E 01, 0.89728244E 01/ 

DATA E/ 0.21 103107E 02, 

2-0 .379597 87E 03 ,-0 .97489220E 02 , 0.12900672E 03, 0. 17949 226E 02, 

2-0 • 129 109 31E 03,-0.557055 746 03, 0. 13524801E 02, 0.14696721E 03, 

2 0.1 79495 2 8E 02 ,-0 . 32 98 1 014E 00, 0.31028836E 02, 0.81657657E 01, 

2 0. 22236961E 02, 0.39124892E 02, 0.81636799E 01/ 

Z=V1 

DO 100 J IM = 1,2 
X = RE A L ( Z ) 

Y=AIMAG(Z) 

E 15 = ( .0,.0) 

AB= CA B S ( Z ) 

IF ( AB.EQ.O. )G0 TO 90 

IF (X.GE.O* .AND. AB.GT.10.)GO TO 80 

YA = ABS ( Y ) 

IF (X .LE .0. .AND. YA.GT.10.JG0 TO 80 

IF IYA-X.GE .17.5.0R.YA.GE .6.5.0R.X+YA.GE.5 .5 .0R.X.GE.3. )G0 TO 20 
IF (X.LE .-9 . )G0 TO 40 
IF IYA-X.GE .2 .5 )G0 TO 50 
IF (X+YA.GE .1.5 ) GO TO 30 
10 N=6.+3.*AB 

E15=1./(N-1.)-Z/N**2 
15 N=N-1 

E 15=1. / (N- 1 . ) -Z *E 15 /N 
IF (N.G6 .3 ) GO TO 15 

E15=Z*E15-CMPLX( .5 772 1 6+ALOG ( AB ) ,ATAN2(Y,X>) 

GO TO 90 
20 Jl=l 
J2 = 6 

GO TO 31 

30 J 1 =7 
J2 = 21 

31 S=(.0,.0) 

YS = Y *Y 

DO 32 I = J 1 , J2 
xi=v U )+X 

CF=W ( I >/<XI*XI+YS> 

32 S=$+CMPLX(XI *CF ,-YA*CF ) 

GO TO 54 

40 T3=X*X-Y*Y 
T4=2.*X*YA 
T5 = X ^T 3- YA *T4 
T6=X£T4+YA*T3 
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Fig. 14a. Subroutine EXPJ 
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UC=CMPtX(D< 11 )+D (12 )*X+0(13)*T3+T5-E( 12)*YA-E< 13)*T4 , 
2 E ( 1 1 )+E ( 12 )*X+E ( 13 )*T3 + T640 ( 12)*YA+D{ 13)*T4> 

V C = CM PLX (D ( 14 J +0 { 15 ) *X+D ( 1 6 )*T3+T5~E ( 15 ) *Y A-E ( 16 )*T4 ? 
2 E ( 14 ) +E ( 15 ) *X + E ( 16 )*T3+T6+0 ( 15 )*YA+D( 16)*T4) 

GO TO 52 
50 T 3=X*X-Y*Y 
T4 = 2 • *X*Y A 


52 

54 


T5=X*T3-YA*T4 
T6=X*T4+YA*T3 
T7=X*T5-YA*T 6 
T8=X*T6+YA*T5 
T9 = X *T 7 -Y A *T 8 


T10=X*T8+YA*T7 

UC = CMPLX(D(1 ) +D { 2 ) *X +D ( 3 ) *T 3+D {4 } **T5 + D(5 )*T7+T9-{ E( 2)*YA+E( 3)*T4 
2 + E (4 ) *T6+E {5 ) *T8 1 *E ( 1 ) +E ( 2 ) *X+E ( 3 ) *T3+E (4)*T5 + Et 5) *T7+T10+ 

3(0(2) *YA + D ( 3 )*T4 +D (4 ) *T6+D (5 )*T8 ) ) 

VC = CMPIX{D(6)+D(7)*X+D(8)#T3+D (9 ) * T5+D ( 10)*T7+T9-(E(7)*YA+E(8)*T4 
2+E (9 )*T6+E (10)*T8) *E ( 6) +E (7 j*X+E <8 )*T3+E (9 )*T5+E< 10 >*T7+T10 + 


3(D(7)*YA+D(8)*T4+D(9)*T6+D( 10 ) *T8 ) ) 
EOUC/VC 
S=EC/CMPLX(X,YA) 

E X = E XP< -X } 

T = EX*CMPLX (COS( YA ) f -SIN (YA) ) 


E 15 =S*T 

5 6 IF (Y.LT.O. )E 1 5 =CDN J G ( E 15) 

GO TO 90 

80 El5=.409319/{ Z + • 19 3 044 )+.421831/(Z+1.02666)+*147126/(Z+2. 56788) + 

2.20633 5E-1/ (Z+4 ,90035 ) + . 1 01 <+ 0 IE -2/ ( Z +8 . 18 2 1 5 ) + . 1 58654E-4/C Z + 
312.7342) +.3 1703 IE -7/ (Z +19 .3957 ) 

E15=E15*CEXP(-Z ) 

90 IF UIM.EG.I )W 1 2 =E 15 
100 Z=V 2 

Z=V2/V 1 

TH = ATAN2( AIMAG(Z ) , RE A L ( Z ) ) - AT AN 2 ( AIM AG < V 2 ) , R E A.L ( V2 ) ) 

2+ ATAN2 ( AIMAGIV1 ) »REAL(V1 ) ) 

A 8 = ASS ( TH ) 

IF ( AB.LT .1 . ) T H = . 0 
IF (TH.GT.1.)TH«6.2831853 
IF (TH.LT .-1*)TH=-6.2831853 
W12=W12-E 15+CMPLX{ •O f TH) 


RETURN 

END 
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Fig. 14b. Subroutine EXPJ 
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The exponential integral Ei(z) is defined in Reference 3. To generate 
W12, subroutine EXPJ calculates Ei(Vl), subtracts Ei(V2) and adds j2mr. 
The term j2n» is determined by the requirement that W12 vanish in the 
limit as VI approaches V2. The integer n may assume values of -1, 0 
or +1. If the integration path does not cross the negative real axis 
in the v plane, n is zero. The term j2n» is calculated below state- 
ment 100. 

APPENDIX 8. Subroutine GANT1 

Subroutine GANT1, listed in Fig. 15, considers the wire structure 
as an antenna. In the input data, VG(J) is the voltage of a generator 
at point IA( J) of segment J. VG(JJ) is the voltage of a generator at 
point IB ( J ) of segment J. The DO LOOP ending with statement 50 uses 
the delta-gap model to determine the excitation voltages CJ(I) for all 
the dipole modes. These are also stored temporarily in CG(I). Then 
subroutine SQROT is called to obtain a solution of the simultaneous 
linear equations. SQROT stores the solution (the loop currents) in 
CJ(I). 

In the DO LOOP ending at statement 80, the complex power input is 
calculated and stored in Yll. GG denotes the time-average power input 
and is the real part of Yll. If the antenna has only one voltage 
generator (with unit voltage and zero phase angle), then Yll also 
denotes the antenna admittance and Zll is the antenna impedance at that 
port. 

Subroutine RITE is called to make the transformation from the loop 
currents CJ(I) to the branch currents CG(J). If IWR is a positive 
integer, RITE will write out the list of branch currents. 

Finally, GANT1 calculates the radiation efficiency EFF. PIN 
denotes the time-average power input. Subroutine GDISS is called to 
obtain the time-average power dissipated. DISS is the total power dis- 
sipated in the lumped loads and the imperfectly-conducting wire. PRAD 
is the time-average power radiated, defined by the difference between 
PIN and DISS. If the antenna has perfect conductivity and purely 
reactive loads, the radiation efficiency is considered to be 100 per 
cent. 

APPENDIX 9. Subroutine SQROT 

Subroutine SQROT is listed in Fig. 16. This subroutine considers 
the matrix equation ZI = V which represents a system of simultaneous 
linear equations. If the square matrix Z is symmetric, SQROT is useful 
for obtaining the solution I with V given. NEQ denotes the number of 
simultaneous equations and the size of the matrix Z. 

On entry to SQROT, S is the excitation column V. On exit, the 
solution I is stored in S. Let Z(I,J) denote the symmetric square 
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SUBROUTINE GANT1 ( I A , IB , INM ,IW R , 1 1 , I 2, I 3, I 1 2 , JA , JB , MO » N, ND »NM , AH OOOl 

2r C » CJ »CG t CMM fO tEFF ,GAM , GG ,CGD , SGD ,VG,Y11,ZU,ZLD,ZS) 0002 

COMPLEX C(1 ) tCJd ) ,CGD«1) tSGO(l ) ,VG(1),ZL0( 1 ) , Y 1 1 1 Z 1 1 , Z S, GAM ,CG( 1) 000 3 

DIMENSION Dd),IAd>,IMl),JA(l),JB(l) 0004 

DIMENSION dC 1) , 12 d ), 1 3(1) ,MD< INM t 4 ) ,NDI 1) 000 5 

2 FQRMATt 1X,1I5»8F10.2) 0006 

5 FO RM AT C 1 HO ) 0007 

DO 50 1*1, N 0008 

CJ(I)=(.0».0) 0009 

K = JA( I ) 0010 

DO 40 KK=1,2 0011 

KA= I A { K ) 0012 

KB= l 8 ( K ) 0013 

J J = K 0014 

F I = 1 * 0015 

IF (KB.EQ.I2d ) )G0 TO 36 0016 

IF (K6.EQ.I1 (I ) ) F I » - 1 . 0017 

C J ( I > =C J ( I ) +F I *V G ( J J ) 0018 

GO TO 40 0019 

36 IF (KA.EQ.I3d ) IF I *-l • 0020 

J J=K+NM 0021 

CJd )=CJ( I )+FI*VG{ JJ) 0022 

40 K=JB(I) 0023 

50 CONTINUE 0024 

DO 55 1*1, N 0025 

55 CG ( I ) =C J ( I > 0026 

CALL SQROT(C,CJtOf I 12, N) 0027 

112=2 0028 

Yll=( .Of .0 ) 0029 

DO 80 1*1, N 0030 

80 YU = Y11 + CJU )*CONJG(CG(I ) ) 0031 

CALL RITE ( IA,IB,INM,IWR,I1, 12,I3tMD,ND,NM,CJ,CG) 00 32 

GG = REAUY11 ) 0033 

Z11 = 1./YU 0034 

P IN =GG 0035 

CALL GD I SS ( AM tCGf CMM t D ,D I SS t GAM ,N M, SGD, ZLD, Z $ ) 00 36 

PRAD= PIN-D I SS 0037 

EFF=100.4PRAD/PIN 0038 

RETURN 0039 

END 0040 

Fig. 15. Subroutine GANT1 
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SUBROUTINE SQROT ( C , S , IW R , 1 1 2 ,NEQ ) OOOl 

complex cm»sm»ss 0002 

2 FORMAT) 1X,1I5,1F10.3»1F15.7,1F 10.0»2F15.6) 0003 

3 FORMAT) 1H0) 0004 

N=NEO 0005 

IF ) 112. EQ. 2 ) GO TO 20 0006 

C(1) = CS0RT(C(1 )) 0007 

00 4 K=2,N 0008 

4 C)K)=C(K)/C(1) 0009 

00 10 I=2»N 0010 

IMO = 1-1 0011 

IP0=I+1 0012 

ID = )I-l)*N-U*I-I)/2 0013 

1 I = ID + 1 0014 

DO 5 1 = 1* IMO 0015 

LI=(L-l)*N-(L*L-L)/2+I 0016 

5 c<m=cm j-c(li)*c(li> 0017 

C( I I )=CSORT(C( II ) ) 0018 

IF ( IPO .GT .N )G0 TO 10 0019 

DO 8 J = I PO * N 0020 

I J = ID + J 0021 

00 6 M= 1 » IMO 0022 

MD=(M-l)*N-(M*M-M)/2 0023 

M I =MD + 1 0024 

M J =MD + J 0025 

6 C 1 I J ) =C ( I J )-C ( M J ) *C ( M I ) 0026 

8 Cl IJ)=C( IJ)/C( II ) 0027 

10 CONTINUE 0028 

20 S ( 1 ) = S ( 1)/C) 1) 0029 

DO 30 1 = 2, N 0030 

IMO = 1-1 0031 

DO 25 1*1 , IMO 003? 

LI=(L-l)*N-(l*L-L)/2+I 0033 

25 S( I )*S( I ) — C C L I )*S(L) 0034 

11=11-1) *N-(I*I-I)/2+I 0035 

30 S(I)*S)I)/Cm) 0036 

NN=( (N + l ) *N ) /2 0037 

S IN ) * S( N )/C (NN ) 0038 

NMO =N- 1 0039 

DO 40 1 = 1, NMO 0040 

K=N- I 0041 

KPU=K+1 0042 

KD=(K-1) *N-(K*K-K)/2 0043 

DO 35 L = KPO ,N 0044 

KL=KD+L 0045 

35 S(K)=S(K)-C(KL)*S(L) 0046 

KK=KD+K 0047 

40 S(K)=S)K)/CCKK) 0048 

IF1IWR.LE.0) GO TO 100 0049 

CNOR= .0 0050 

DO 50 1=1, N 0051 

SA = CABS ( S ( I ) ) 0052 

50 IF (SA.GT.CN OR )CNOR=SA 0053 

IF (CNOR.LE .0. )CN0R=1 . 0054 

DO 60 1 = 1, N 0055 

SS=S(I) 0056 

SA = CAB S ( SS ) 0057 

SNOR= SA/CNO R 0058 

PH=.0 0059 

IF ( S A .GT .0 . )PH=57.29578*ATAN2( AIMAG(SS) ,REAL(SS) ) 0060 

60 WRITE (6,2 )I ,SNOR,SA, PH, SS 0061 

WRITE 16,3) 006? 

100 RETURN 00 63 

END 006A 

Fig. 16. Subroutine SQROT 
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matrix. On entry to SQROT, the upper-right triangular portion of Z(I,J) 
is stored by rows in C(K) with 

(4) K = (I - 1 )*NEQ - (1*1 - I) / 2 + J 

If 112 = 1, SQROT will transform the symmetric matrix into the auxiliary 
matrix (implicit inverse), store the result in C(K) and use the aux- 
iliary matrix to solve the simultaneous equations. If 112 = 2, this 
indicates that C(K) already contains the auxiliary matrix. 

The transformation from the symmetric matrix to the auxiliary matrix 
is programmed above statement 10, and the solution of the simultaneous 
equations is programmed in statements 20 to 40. If IWR is positive, the 
program below statement 40 will write the solution. 

SQROT uses the square root method described in Reference 4. The 
original symmetric matrix Z and the upper triangular auxiliary matrix A 
are related by 

(5) Z = A' A 

where A' is the transpose of A. 

In the thin-wire application, SQROT must be called with 112 = 1 
before it is called with 112 = 2. With a large matrix, the execution 
time in SQROT is much smaller with 112 = 2 than with 112 = 1. 

APPENDIX 10. Subroutine RITE 

Subroutine RITE is listed in Fig. 17. Given the list of loop 
currents CJ(I), this subroutine generates a list of branch currents 
CG(J) . CG(J) and CG(JJ) denote the currents at IA(J) and IB ( J ) , re- 
spectively, on the wire segment J, where JJ = J + NM ; If IWR is a 
positive integer, the program below statement 110 writes a list of the 
branch currents. The symbols in this list are defined as follows: 


K 

the segment number 

at IA(K) 

ACJ 

normalized current magnitude 

BCJ 

normalized current magnitude 

at I B ( K) 

PA 

phase of current at IA(K) 


PB 

phase of current at IB( K) 


CJA 

complex current at IA(K) 


CJB 

complex current at IB(K) 



The phase angles PA and PB are in degrees. Even if IWR is negative, 
RITE generates the branch-current list for use in subroutine GDISS. 
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SUB RUU TINE R I TE ( I A , I B , I NM , IW R , I 1 , I 2, I 3, MD, ND, NM ,C J , CG ) OOOl 

COMPLEX CJtl>,CG(l),CJA,CJB 0002 

DIMENSION IA(1),IB<1), 11(1), 12(1), 13(1), MDIINM, 4), ND(1> 0003 

2 FORMAT ( IX, 1 15 »2F10.3,2F10.0»4F15»6) 0004 

5 FORM AT 1 1 HO ) 0005 

AM AX = . 0 0006 

DO 100 K= 1 ,NM 0007 

KA= I A ( K) 0008 

KB= I B ( K ) 0009 

CJA=(.0,.0) 0010 

CJB=(.0,.0) 0011 

NDK=ND ( K ) 0012 

DO 40 11=1, NDK 0013 

I=MD!K,II) 0014 

F I = 1 . 0015 

IF (KB.E0.I2 ( I ) JGO TO 36 0016 

IF (KB.EQ.Iim )FI«-1. 0017 

CJA=CJA+F J+CJII) 0018 

GO TO 40 0019 

36 IF (KA.EQ.I3U) IF I=-l. 0020 

CJB=CJB+F I*CJ(I) 0021 

40 CONTINUE 0022 

CG ( K ) =C JA 0023 

KK=K+NM 0024 

CG ( KK ) = CJ B 0025 

AC J=C A8S ( C J A ) 0026 

BC J=CABS ( CJB ) 0027 

IF <ACJ.GT.AMAX)AMAX=ACJ 0028 

IF (BCJ.GT.AMAX)AMAX=BCJ 0029 

100 CONTINUE 0030 

IF ( iy R .GT .0 JGO TO 110 0031 

RETURN 0032 

110 IF (AMAX.LE .0. JAMAX=1. 0033 

DO 200 K=1,NM 0034 

C JA=CG( K ) 0035 

KK=K+NM 0036 

CJB=CG(KK) 0037 

ACJ=CABS(CJA)/AMAX 0038 

BCJ=CABS(CJB)/AMAX 0039 

PA=57.29578*ATAN2 (AIMAG(CJA) ,RE ALICJA) ) 0040 

PB = 5 7.295 78*ATAN2(AIMAG(CJB) , REAL I CJB) ) 0041 

200 WRITE (6,2)K,ACJ,BCJ,PA,PB,CJA,CJB 0042 

W R I TE ( 6 ,5 ) 0043 

RETURN 0044 

END 0045 


Fig. 17. Subroutine RITE 
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APPENDIX 11. Subroutine GDISS 

Subroutine GDISS is listed in Fig. 18. This subroutine uses 
Eq. 50 of Reference 1 to calculate the time-average. power. dissipated 
in the imperfectly conducting wire. This is accomplished in the DO 
LOOP terminating at statement 100. The power dissipated in the lumped 
loads is calculated in the DO LOOP terminating with statement 140. 

DISS denotes the time-average power dissipated in the wire and the 
1 oads . 

APPENDIX 12. Subroutine GNFLD 

Subroutine GNFLD, listed in Fig. 19, inputs the loop currents 
C J ( I ) , calls GNF for the near-zone field of each wire segment, and 
sums over all the segments to obtain the near-zone field of the wire 
antenna or the near-zone scattered field of the wire scatterer. EX, 

EY and EZ denote the cartesian components of this field at the obser- 
vation point (XP.YP.ZP). This calculated field does not include the 
incident fields of the magnetic frills or loops associated with gener- 
ators on the antenna. It also does not include the radiation from the 
polarization currents in the dielectric insulation. 

This subroutine could be simplified and speeded by inputting the 
branch currents CG(J) instead of the loop currents CJ(I). However, 
this would increase the storage requirements because the far-field 
subroutine GFFLD would have to store the branch currents induced by 
the phi -pol ari zed and theta-polarized incident waves. 

APPENDIX 13. Subroutine GNF 

Subroutine GNF, listed in Fig. 20, uses Eqs. 75 and 76 of . 
Reference 1 to calculate the near-zone electric field of a sinusoidal 
electric monopole with endpoints at (XA,YA,ZA) and (XB,YB,ZB). The 
observation point is at (X,Y,Z). EX1, EY1 and EZ1 are the components 
of the field generated by the mode with unit current at (XA,YA,ZA). 

EX2, EY2 and EZ2 denote the field generated by the mode with unit 
current at (XB,YB,ZB). GNF is similar to GGS, and Appendix 5 defines 
many of the symbols used in both subroutines. 

APPENDIX 14. Subroutine GFFLD 

The far-field subroutine GFFLD, listed in Fig. 21, is discussed in 
section II. In antenna gain calculations with INC = 0, the loop, 
currents CJ(I) are employed by GFFLD to calculate the far-zone field. 
The field of each segment is obtained by calling GFF, and a surranation 
over all the segments yields the field of the antenna. 

In a bistatic scattering situation with INC = 2, the input data 
include the loop currents EP and ET induced by phi -pol ari zed and theta- 
polarized incident waves. These currents were calculated by GFFLD in a 
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SUBROUTINE GO I SS ( AM f CG , CMM ,D ,DISS,GAM ,N M f SGD f ZLD,Z S) 
COMPLEX CG(1 ),SGOm f ZLO(l J »CJA f CJB,GAM,ZS 
DIMENSION 0(1) 

DATA PI/3* 14 159/ 

D I SS= • 0 

IF (CMM • LE .0. )G0 TO 120 
ALPH= RE A L ( GAM ) 

BE T A = A IM AG ( GAM ) 

RH=REAL(ZS)/(4.*PI*AM) 

DO 100 K= 1 ?NM 
DK=D ( K ) 

DEN=CABS(SGO(K) )**2 
EAD=EXP( ALPH*DK) 

CAD=(E AD+l./EAD )/2. 

CBD = COS ( BE T A*DK ) 

SAD =DK 

IF (ALPH.NE *0*)SA0=(EAD-1./EAD)/(2**ALPH) 

SBD =DK 

IF ( BE TA *NE .0 . ) SBO =S IN ( BE T A*D K > / BE T A 
F A = RH*( SAD*CAD-SBD*CBD J/DEN 
FB = 2 **RH*(CAD*SBD-SAD*CBD ) /DEN 
CJ A=CG ( K ) 

L=K+NM 
C J B = CG ( L) 

100 D ISS=DI SS+F A*(CABS(CJA ) **2 +C ABS ( CJ B )**2 ) 

2+FB* (REAL(CJA)*REAL(CJB )+AIMAG(C JA >* AI MAG ( C JB ) ) 

120 DO 140 J- 1 »NM 
K = J +NM 

14 0 D ISS=DISS+REAL(ZLD( J > )*(CABS(CG( J) )**2) 

2+REAUZLD (K) )* (CABS(CG( K ) )#*2 ) 

RETURN 

END 
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Fig* 18. Subroutine GDISS 
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SUBROUTINE GNF LD ( I A , I B , INM , 1 1 , 1 2 , 1 3 >M 0,N ,ND,NM, AM, CGD, SGO, ETA, GAM 
2»CJ*D,X,Y,Z,XP,YP,ZP,EX,EY,EZ) 

COMPLEX E X ,EY,EZ ,EX1 ,EY1 ,EZ1,EX2,EY2 ,EZ2,ETA,GAM 
. COMPLEX Cjm,CGD(l>,SGO<n 

DIMENSION IA ( 1 ) , IB(1), 11(1), 12(1), 13(1), Dll), X(l), Y(l), Z ( 1) 
DIMENSION MDIINM ,A ) ,ND Cl) 

DATA PI , TP/3. H159 ,6.28318/ 

EX= ( .0, .0) 

EY=( .0,.0) 

EZ=C.0,.0) 

DO 140 K= 1 «NM 
KA= I A ( K ) 

K8= I B ( K ) 

CALI GNMX(KA) »Y(KA) tZ(KA) f X { KB > t Y < KB ) ? Z ( KB ) tXPfYP*ZP*AM f D(K) 

2 » C GD { K ) ,SG0(K)tETA f GAM,EXl,EYl f EZl»EX2»E Y2tEZ2) 

NDK=NDU) 

DO 140 1 1 = 1 1 ND K 
I=MD ( K , 1 I ) 

F 1 = 1 . 

IF ( KB *E0 • 12 ( I ) ) GO TO 136 
IF<KB.E0.I1(I))FI»-U 
E X=E X+F I *E X 1 *C J ( I ) 

EY=EY+FI*EY1*CJ(I ) 

EZ=EZ+F I*EZ1*CJ< I ) 

GO TO 140 

136 IFtKA.EQ.I3U ) )F I =-l • 

EX=EX+FI*EX2*CJ(I ) 

E Y=E Y+F I *E Y2*CJ ( I ) 

E Z=E Z + F I +E 22 *CJ ( I ) 

140 CONTINUE 
RETURN 
END 
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Fig. 19. Subroutine GNFLD 
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.SUMII/Iil JM- GNF (XA,YA,ZA,XB,YB,ZB,X,Y ,Z , AM, DS, CGDS, SGDS , E TA, GAM 0001 

2,EX1,LY1,EZ1,EX2,EY2,EZ2> 0002 

COMPLEX F JA,FJB,F J1,E J2.ER1 ,ER2 ,E S1,E S2, SGDS, GAM, CS T.CGOS, ETA 0003 

COMPLEX EX1 ,EY1,EZ1 ,EX2,EY2,EZ2 0004 

DATA PI/3.14159/ 0005 

CAS= ( XB-XA )/DS 0006 

CBS= ( YB-YA )/DS 0007 

CGS= ( ZB-ZA ) /D S 000B 

SZ=(X-XA)*CAS+(Y-YA)*CBS+(Z-ZA)*CGS 0009 

ZZ1=SZ 0010 

ZZ2=SZ-DS 0011 

XXZ=X-XA-SZ*CAS 0012 

YYZ=Y-YA-SZ*CBS 0013 

ZZZ=Z-ZA-SZ*CGS 0014 

RS=XXZ**2+YYZ**2+ZZZ**2 0015 

R1=SQRT (RS+ZZ1**2) 0016 

EJA=CEXP(-GAM*R1) 0017 

EJ1=EJA/R1 0018 

R2=SQRT(RS+ZZ2**2> 0019 

EJB=CEXP(-GAM*R2) 0020 

E J2=E J B/R2 0021 

ES1=EJ2-EJ1*CGDS 0022 

ES2=EJ1-E J2*CGDS 0023 

ER1=(.0,.0) ' 0024 

ER2=( . 0 , . 0 ) 0025 

AM S = AM *AM 0026 

IF ( RS. LT. AMS )G0 TO 80 0027 

CTH1=ZZ 1/R1 0028 

CTH2=ZZ2/R 2 0029 

E R 1= ( E JA*SGOS+E JA#CGDS*CTH1-E JB*CTH2)/RS 0030 

ER2=(-E JB*SGDS+E JB*CGDS*CTH2-E JA*CTH1)/RS 0031 

80 CST=ETA/(4.*P1*SGDS) 0032 

EX1=CST*(ES1*CAS+ER1*XXZ > 0033 

EY1=CST*(ES1*CBS+ER1*YYZ) 0034 

EZ1=CST*(ES1*CGS+ER1*ZZZ ) 0035 

EX2=CST*(ES2*CAS+ER2*XXZ) 0036 

EY2=C5T*(E52*CBS+ER2*YYZ ) 0037 

EZ2=CST*(ES2*CGS+ER2*ZZZ> 0038 

RETURN 0039 

END 0040 


Fig. 20. Subroutine GNF 
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SUBROUTINE GF F LD C I A , IB , INC »INM , IW R , II « I 2 , I 3, 1 1 2, MO ,N , NO .NM, AM OOOl 

2, ACSP,ACST,C,CGD,CG,CJ,CMM,D,ECSP,ECST,EP,ET,EPP,ETT,EPPS,EPTS 0002 

3, ETPS,ETTS *GG,GPP,GTT ,PH,SGD , SCSP , SCS T, SPPM , SP TM, STPM, S TTM, TH 0003 

4, X,Y,Z,ZLQ,ZS,ETA,GAM) 0004 

COMPLEX CJI ,ET1,ET2,EP1,EP2,EPPS,ETTS,EPTS,ETPS,ZS,VP,VT 0005 

COMPLEX Cm,Cjm,EP(l),ET(l),ePP(l>,ETT(l),ZLOm 0006 

COMPLEX ETA,GAM,CG0(1) ,SGO ( 1 J ,CG< 1 ) 0007 

DIMENSION I A (1) , 1 B(l) ,1 1< 1) ,I2( 1) ,I3( 1 ) »ND( 1) ,M0< INM,4) 0008 

DIMENSION D(1),X(1),Y(1),Z(1) 0009 

DATA PI , TP/3. 14159 ,6.28318/ 0010 

CJI=-4.*PI/(ETA*GAM) 0011 

GGG = REAL( 1 ./ETA) 0012 

THR=.0174533*TH 0013 

CTH=COSCTHR) 0014 

STH=SIN I THR ) 0015 

PHR =. 0174533 *PH 0016 

CPH=COS(PHR) 0017 

SPH= S IN ( P HR ) 0018 

DO 130 1 = 1, N 0019 

ETTI 1) = ( . 0 , .0 ) 0020 

130 EPP(I )=( .0, .0) 0021 

DO 140 K= 1 ,NM 0022 

KA= I A IK) 0023 

KB= 1 B I K ) 0024 

CALL GFFIX(KA) ,Y(KA) , Z ( K A ) ,X(KB),Y(KB),Z(KB) ,D(K) 0025 

2,CGD<K) ,SGO«K) ,CTH ,STH, CPH.SPH ,G AM ,E TA,E T1,ET2,EP1,EP2) 0026 

NDK=ND( K ) 0027 

DO 140 11=1 ,NDK 0028 

I =MD ( K , 1 1 ) 0029 

F I =1 . 0030 

IF (KB.EQ.I2 ( I ) )G0 TO 136 0031 

IF IKB.EQ.IIU) )FI=-1. 0032 

EPPl I )=EPP( I ) +F I*EP1 0033 

ETTI I >=£TT(I )+F I *ET1 0034 

GO TO 14 0 0035 

136 IF ( KA.EQ . 131 I ) IF I =-l . 0036 

EPP1I )=EPP(I)+FI*EP2 0037 

ETT ( I ) =ETT ( I ) +F I *ET2 0038 

140 CONTINUE 0039 

E PPS= ( .0 , .0 ) 0040 

ETTS= ( .0, .0 ) 00^1 

IF I INC .EQ .0 )G0 TO 200 0042 

IF i INC.EQ.2 )G0 TO 170 0043 

DO 150 1 = 1, N 0044 

ET I I )=ETT ( I )*CJ I 0045 

150 E P< r >=E PP ( I ) *CJ I 0046 

CALL SQR0T(C,EP»0,112,N) 0047 

112=2 0048 

CALL SQROT (C,ET ,0,112,N ) 0049 

CALL RITE <IA, IB, INM.IWR, 11,12, I3,MD,ND,NM,EP,CG) 0050 

CALL GDI SS< AM ,CG,CMM ,D ,P0 IS, GAM ,NM , SGD,ZLD,ZS) 0051 

CALL RITE ( I A ,1 B » INM , IWR, II, 12,13, MD,ND,NM,ET,CG) 0052 

CALL GO 1 SS ( AM ,CG,CMM ,D ,TD I S ,GAM ,NM ,SGD,ZLD,ZS) 0053 

ACSP=PD IS/GGG 0054 

ACST=TD I S/GGG 0055 

PIN= .0 0056 

T IN= .0 0057 

DO 164 I = 1 *N 0058 

VP=CJI*EPP(I) 0059 

VT = CJ I *E TT < I ) 0060 

P IN = PIN + RE AL ( V P*CON JGI E P 1 1 ) > ) 0061 

164 TIN=T IN + REAL (VT *CONJG(ET ( I ) ) ) 0062 

Fig. 21a. Subroutine GFFLD 
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ECSP-PIN/GGG 0063 

ECST*T IN/GGG 0064 

SCSP=ECSP-ACSP 0065 

SCST=ECST-ACST 0066 

170 E PTS= ( • 0 , .0 ) 0067 

ETPS=<.0,.0) 0068 

DO 180 1 = 1, N 0069 

EPPS=EPPS+EP(I )*EPP(I ) 0070 

EPTS=EPTS+EP( I )*ETT( I) 0071 

ETTS=ETTS+ET(I >*ETT< I ) 0072 

180 ETPS=ETPS+ET( I )*EPP( I > 0073 

SPPM=2.*TP*(CABS<EPPS)**2) 0074 

SPTM=2.*TP*(CA8S(EPTS) ** 2 ) 0075 

STPM=2.*TP*(CABS(ETPS)**2) 0076 

STTM=2.*TP*(CABS(ETTS)**2 > 0077 

RETURN 0078 

200 0° 260 1 = 1, N 0079 

E TT S=E TTS + CJ 1 1 ) *E TT ( I ) 0080 

260 EPPS=EPPS+CJ( 1 )*EPP( 1 ) 0081 

APP=CABS(EPPS) 0082 

ATT = CAB S ( E TTS ) 0083 

GPP=4.*PI *APP*APP*GGG/GG 0084 

GTT=4 .*P I *ATT *ATT *GGG/GG 0085 

RETURN 0086 

END 0087 


Fig. 21b. Subroutine GFFLD 
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previous call for the backscattering situation with INC = 1. Thus, 
a bistatic call must be preceded by a backscatter call. 

EPP(I) and ETT(I) denote the phi-polarized and theta-polarized 
far-zone fields of dipole mode I with unit terminal current. In a 
backscattering situation, the excitation voltages EP(I) and ET(I) are 
obtained by multiplying EPP and ETT by the constant CJI. (See Eqs. 

38 , 39 and 40 in Reference 1.) Then calls are made to SQROT which 
stores the solution (the induced loop currents) in EP(I) and ET(I). 

RITE is called for the branch currents CG(J), and GDISS is called for 
the time-average power dissipated in the imperfectly conducting wire 
and the lumped loads. This power is denoted PDIS and TDIS for phi - 
polarized and theta-polarized incident waves, respectively. 

In scattering problems, the incident plane wave has unit electric 
field intensity at the coordinate origin. GGG denotes the time-average 
power density of the incident wave at the origin. ACSP and ACST denote 
the absorption cross sections for the phi and theta polarizations. 

PIN and TIN denote the time-average power input to the wire struc - 
ture, delivered by the equivalent voltage generators VP and VT at the 
terminals. PIN and TIN apply for the phi and theta polarizations, 
respectively. The time-average power input is regarded as the sum of 
the time-average power dissipated (in the wire and the lumped loads) 
and the time-average power radiated or scattered by the wire. ECSP 
and ECST denote the extinction cross sections and SCSP and SCST are 
the scattering cross sections. 

The distant field is calculated in the DO LOOP ending with state- 
ment 180 for scattering situations, ?nd in the DO LOOP ending with 
statement 260 for the antenna situation. In these fields, the range 
dependence is suppressed as in Eq. (1). 

The radar cross sections (echo areas) SPPM, SPTM, STPM and STTM 
are defined as in Eq. 72 of Reference 1 with the incident power density 
(S,- or GGG) evaluated at the coordinate origin. The user selects the 
location of the origin when supplying the input data for the coordinates 
of all the points on the wire. 

For an antenna, the following definition is employed for the power 

gain: 

(6) Gp(8,4>) = limit 4irr 2 e 2ar S(r,e,4>) / P n - 

where P^ (or GG in the program) denotes the time-average power input 
and S(r,9,<f>) is the time-average power density in the radiated field. 

For an antenna in a lossless medium, a vanishes and Eq. (6) reduces to 
the standard definition of power gain. Without the factor e 2otr in 
Eq. (6), the power gain would vanish for a finite antenna in a conduct- 
ing medium. GPP and GTT denote the power gains associated with the phi - 
polarized and theta-polarized components of the field, respectively. 
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APPENDIX 15. Subroutine GFF 

Subroutine GFF, listed in Fig. 22, uses the equations in Appendix 2 
of Reference 1 to calculate the far- zone field of a sinusoidal electric 
monopole. The monopole has endpoints (XA,YA,ZA) and (XB,YB,ZB). 

EP1 and ET1 denote E. and E^ for the mode with unit current at (XA,YA,ZA). 
EP2 and ET2 denote the fields for the mode with unit current at 
(XB,YB,ZB). The range dependence is suppressed as in Eq. (1). The 
far field vanishes in the endfire direction where GK = 0. 


49 



SUBROUTINE GF F ( X A , Y A , Z A , XB , Y B , Z B ,D , 0001 

2CGO.SGD ,CTH,STH,CPH, SPH, 0002 

2GAM .ETA.ETl ,E T2 ,EP1 ,EP2 ) 0003 

COMPLEX E T1 ,ET2 ,E PI ,E P2 ,GAM ,ET A 0004 

COMPLEX GO tCGD » SGD ?E GO 0005 

COMPLEX EGFA,EGFB,EGGD,ESA,ESB 0006 

COMPLEX CST 0007 

FP=12. 56637 0008 

XAB= XB-XA 0009 

YAB=YB-YA 0010 

Z AB= ZB-Z A 0011 

CA = XAB/D 0012 

CB=YAB/D 0013 

CG=Z A3/D 0014 

G=(CA*CPH+CB*SPH)*STH+CG*CTH 0015 

GK=1.-G*G 0016 

ET1=(.0,.0) 0017 

ET2 = < .0 i .0 ) 0018 

EPl=(.Ot.O) 0019 

EP2=(.0,.0) 0020 

IFIGK.LT. .OODGO TO 200 0021 

F A=(XA*CPH+YA*SPH )*STH+Z A*CTH 0022 

FB= (XB*CPH+YB*SPH)*STH+ZB*CTH 0023 

EGFA = CEXP( GAM *F A ) 0024 

EGFB = CEXP(GAM*FB) 0025 

EGGD=CE XP (GAM*G*D ) 0026 

CST=ETA/<GK*SGD*FP) 0027 

ESA=CST*E GF A * [ E GGD-G*SGD -CGD ) 0028 

ESH = CST--EGFB*( 1 . /E GGD +G*S GO -CGO ) 0029 

T=(CA*CPH+CB 4 SPH)*CTH-CG*STH 0030 

P = -C A* SPH+CB *C PH 0031 

E T1 = T*E SA ' 0032 

ET2=T#ESB 0033 

E P1 = P*E SA 0034 

EP2=P*ESB 0035 

200 CONTINUE 0036 

RETURN 0037 

END 0038 


Fig. 22. Subroutine GFF 
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